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0. Abstract 



Bounds on cross section measurements of chargino pair production at LEP yield a boimd on the 
chargino mass. If the GUT relation Mi =5/3 tan^ B-a)M2 is assumed, then the lightest neutralino 
must be heavier than 45 — 50 GeV. If Mi is considered as a free parameter independent of M2 
there is no bound on the mass of the lightest neutralino. In this thesis, I examine consequences 
of light, even massless neutralinos in cosmology and particle physics. 

In Chapter |2l I discuss mass bounds on the lightest neutralino from relic density measure- 
ments. The relic density can be calculated by solving the Boltzmann equation. If the relic den- 
sity is considered as a function of the particle mass then there are two mass regions where 
the relic density takes on realistic values. In the first region the neutralino is relativistic and 
its mass must be lower than 0.7 eV, in the second region the neutralino is nonrelativistic and 
must be heavier than 13 GeV. I compare the Cowsik-McClelland bound, the approximate 
solution of a relativistic particle for the Boltzmann equation, and the Lee- Weinberg bound, the 
non-relativistic approximation, with the full solution and I find that the approximation and the 
full solution agree quite well. 

In Chapter |3l I derive bounds on the selectron mass from the observed limits on the cross sec- 
tion of the reaction e+e^ X.1X.2 LEP, if the lightest neutralino is massless. If M2, }i < 200 GeV, 
the selectron must be heavier than 350 GeV. 

In Chapter m I study radiative neutralino production e+e^ — > XiXil the linear collider 
with longitudinally polarised beams. I consider the Standard Model background from radia- 
tive neutrino production e+e^ — > vvy, and the supersymmetric radiative production of sneu- 
trinos e+e^ —>■ vv*^, which can be a background for invisible sneutrino decays. I give the 
complete tree-level formulas for the amplitudes and matrix elements squared. In the Minimal 
Supersymmetric Standard Model, I study the dependence of the cross sections on the beam po- 
larisations, on the parameters of the neutralino sector, and on the selectron masses. I show that 
for bino-like neutralinos longitudinal polarised beams enhance the signal and simultaneously 
reduce the background, such that search sensitivity is significantly enhanced. I point out that 
there are parameter regions where radiative neutralino production is the only channel to study 
SUSY particles at the ILC, since heavier neutralinos, charginos and sleptons are too heavy to be 
pair-produced in the first stage of the linear collider with a/s = 500 GeV. 

In Section I4.4[ I focus on three different mSUGRA scenarios in turn at the Higgs strahlung 
threshold, the top pair production threshold, and at a/s = 500 GeV. In these scenarios at the cor- 
responding i/s, radiative neutralino production is the only supersymmetric production mech- 
anism which is kinematically allowed. The heavier neutralinos, and charginos as well as the 
sleptons, squarks and gluinos are too heavy to be pair produced. I calculate the signal cross 
section and also the Standard Model background from radiative neutrino production. For my 
scenarios, I obtain significances larger than 10 and signal to background ratios between 2% 
and 5%, if I have electron beam polarization P^- — 0.0 — 0.8 and positron beam polarization 
Pg+ = 0.0 — 0.3. If I have electron beam polarization of Pg- = 0.9, then the signal is observable 
with Pg+ = 0.0 but both the significance and the signal to background ratio are significantly 
improved for Pg+ = 0.3. 
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0. Abstract 



In Chapter |5l I present a method to determine neutralino couplings to right and left handed 
selectrons and Z bosons from cross section measurements of radiative neutralino production 
and neutralino pair production e^e^ /^i/t2/3/4' ~^ X2X.2 ll^C. The error on the cou- 

plings is of order 0(0.001 — 0.01). From the neutralino couplings the neutralino diagonalisation 
matrix can be calculated. If all neutralino masses are known. Mi, M2, and ji can be calculated 
with an error of the order 0{1 GeV). If also the cross sections of the reactions e+e^ ^2^3/4 
can be measured the error of Mi, M2, and }i reduces to 0{1 GeV). 
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1. The Gaugino sector in the MSSM 



The Standard Model (SM) has been tested to high precision. But many problems remain un- 
solved. The SM does not include gravity. The electro-weak couplings and the strong coupling 
do not unify in a point at the GUT scale Aqut- The SM model cannot explain why there is 
so much more matter than antimatter in the universe and it does not provide a dark matter 
candidate. 

One solution to these problems is supersymmetry [1-5]. In supersymmetric theories, each 
fermion is mapped onto a boson and vice versa. The spin of the fermion and its partner boson 
differ by half a unit, the other quantum numbers are unchanged. 

The superpartners of leptons, quarks, gauge bosons, and Higgs bosons are called sleptons, 
squarks, gauginos, and higgsinos, respectively. The two neutral gauginos Aq, A3 and the two 
neutral higgsinos h\, h\ have the same quantum numbers and mix. The physical mass eigen- 
states are obtained by diagonalisation of the mass matrix. These neutral particles are called 
neutralinos. They are Majorana fermions. The two charged gauginos and two charged higgsi- 
nos mix to charginos. 

At low energies no SUSY particles have been observed, so SUSY must be broken. The most 
common way is introducing explicitly soft SUSY breaking terms. 

The part of the Lagrangian which describes the neutralino mixing is given by [4] 



C 



1 

-AoAoMi 



-A3A3M2 



vihl) 



(1.1) 



with 



M 



( 



Ml 






M2 



— mz sin 0i„cos/3 
\ m-z sin Q-^o sin /S 



mz cos Q-u, cos f> 
-mz cos Qio sin/3 



-?«zsin0i„cos/6 
mz cos Qio cos /3 




mz sin 0j(,sinj6\ 
—mz cos Qio sin/5 

/ 



= I An, A3 



(the are Weyl spinors). 



(1-2) 



(1.3) 



Ml and M2 are the U{\)y and the S!J(2)a, gaugino mass parameters, respectively. They break 
SUSY explicitly. }i is the higgsino mass parameter and tan /6 = ^ is the ratio of the two vacuum 
expectation values of the Higgs fields, mz the Z boson mass, and tan 9^^ the weak mixing angle. 
Ml, M2, and }i are real parameters, if CP is conserved, in general they are complex: 



Ml 



|Mi|e 



(1.4) 



The matrix M is symmetric, even for M complex. The reason for this fact is that in Eq. 11.11 
there appear no hermitian conjugated fields. The matrix M can be diagonalised by an unitary 
matrix N using Takagi's factorization theorem [6] 



diag(m^o,m^o,m^o,m^o) = N*MN 



(1.5) 
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1. The Gaugino sector in the MSSM 



The diagonal elements m^o are non-negative and are the square roots of the eigenvalues of 

The transformation Eq. (|1.5|) is not a similarity transformation, if N is complex. 

If M is a real matrix it can also be diagonalised by an orthogonal matrix. From the lower right 
2x2 submatrix one can see that at least one eigenvalue is negative. This sign is interpreted as 
the CP eigenvalue of the neutralino. The masses of the neutralinos are \mi\,i = 1 ... 4. The sign 
can be absorbed in a phase of the corresponding eigenvector, leading back to Eq. (|1.5|) . 

The eigenvalues of MM+ and the diagonalisation Matrix N can be obtained algebraically, see 
Ref. [7] or numerically. The algebraic method is problematic because it is numerically unsta- 
ble. Gunion and Haber present in [8] approximate solutions to the eigenvalues of M and the 
diagonalisation matrix N, if |Mi^2 i mI ^ "^z- 

Without loss of generality M2 can be chosen positive. Proof: Let M2 = |M2|e"''2. The phase (p2 
of M2 can be removed by the transformations: 

/Ao\ 



A3 



1^' 



/Aoe-^'^/2^ 
h\ ei^2/2 



(1.6) 



The parameters Mi and }i transform as 

Ml 



M[ = Ml e^ 



(1.7) 

^■=^e-'1'^. (1.8) 

It is not necessary to transform the higgsino fields. Alternatively, the vacuum expectation values 
^1/2 can be transformed as — '01/2^^^'^^^'^ leaving tan j6 invariant. 

li (p2 = 71 then the signs of Mi and ji are interchanged. This transformation reverses also the 
sign of the eigenvalues of M. 

In GUT theories. Mi and M2 are related by 



Ml 



-tan^ 



?i«M2 ~ -M2 . 



(1.9) 



It follows that Ml and M2 can be chosen positive 
M can have zero eigenvalues. From det(M) : 
Ref. [9], 

ji [M2m\ sir? e,„ sin(2|6) + Mi ( - M2}i 
M2m| sin^ 0i„ sin(2;6) 



= det(M) 



it follows in the CP conserving case, see 
m|cos2 0,„sin(2jS))] 

^ ~ M2}i - m\ cos2 Q,„ sin(2|6) ' ^^'^^^ 
is excluded due to experimental constraints from the Z^-widths measured 



V M 



The solution ^ 
at LEP [10]. 

In the CP violating case, substitute Mi 1-^ Mie^*^!, }i 
equations, which must be separately fulfilled: 



^e"^'' with Ml, f/ > 0. This yields two 



D^m det(M) =0 => ]i = 
5Kedet(M)=0 ^ Mi 



m| cos^ Qio sin(2/3) sin (\)\ 
M2sin((/)i +^,,) 

sin(^i +^,,) 



-M2 tan^ Qu 



(1.11) 



(1.12) 
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or 



m\ cos^ d-,u sin(2^) sin <pi m\ sir? O-a, sin(2/3) sin (pi 

M2 = : — 7 ^ ana Mi = : , [^-^j) 

or 

. ^M2sm{(pl + (p^) 2fl sin((/)i + <^j,) 

smfzHj = — 7y r — and Mi = — M2tan t^n, ^ —. (1-14) 

?«2 cos^ 9zi, sin (/)i sin 

It follows immediately that sint/^i/ sin^^, < and sin((/)i + (p^) / sirv(p^ < must hold. Also 
in the CP violating case one can always find parameters \Mi\,(pi, M2, \ }i\,(pii, and tan^ to get 
m^o = 0. 

The chargino mixing is described by the following matrix: 

C = -{xp-fXxp+, (1.15) 

/ M2 Vlniw sin B\ 

X = , (1.16) 

\V2mwcos/S ji ) 

xp+ = {A+,hlf, ^p- = {\-,h\Y, (1.17) 

X is not symmetric, so it must be diagonalised by a biunitary transformation: 

diag{m^,m^) = ^J*Xy-^ (1.18) 

with U, V unitary 2x2 matrices. The matrices U and V are obtained by solving 

diag {{m^f, (^2 )^) = VX+XV-^ = U*XX+U^ . (1.19) 

The eigenvalues can be obtained analytically, see Ref. [1,2]. In practical use it is easier to diago- 
nalize the matrix X numerically but using the analytical formulae. 

The lower experimental bound on the lightest chargino mass is m^± > 104 GeV [11]. This 
bound leads to lower bounds on /y and M2: }i,M2 > 100 GeV. If Eq. ( 11.91 ) is assumed, then Mi 
depends on M2 and from this fact follows a lower bound on the mass of the lightest neutralino: 
m-o > 49 GeV [12]. But up to now there is no evidence that Eq. (|1.9|) holds. So I consider Mi 
as a free parameter. In the following I study implications of massless and light neutralinos. In 
Chapter 121 I discuss bounds on the neutralino mass from dark matter density measurements. 
In Chapter |3l I derive bounds on the selectron mass from the observed cross section limits 
from X.1X2 production at LEP, if Xi is massless. In Chapter HI I calculate the cross section for 
radiative neutralino production and its neutrino and sneutrino background at a future e+e^ 
linear collider. I discuss the influence of beam polarisation on radiative neutralino production 
and consequences of SUSY searches at a future linear collider. Finally, in Chapter |5l I present a 
method how to determine neutralino couplings to the right and left handed selectron and the Z 
boson. 
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2. Cosmological bounds on neutralino masses 



2.1. The Cowsik-McClelland-bound 

I derive bounds on the mass of the lightest neutralino through cosmological considerations. 
Neutralinos are neutral and interact only weakly. If they are (pseudo-)stabile, they are dark - 
matter (DM) candidates. The dark matter density Qdm^^ has been measured by the WMAP col- 
laboration [13]. This constrains the mass(es) of the particle(s) which constitute the dark matter. 

In Ref. [14], Kolb and Turner describe the thermal evolution of the Universe and its impact 
on particle physics. I give a short summary in order to clarify the subsequent section. 



2.1.1. The Expansion of the Universe 

The expansion of the Universe is described by the Einstein field equations with the Robertson- 
Walker (RW) metric. In the RW metric, the Universe is assumed to be homogeneous and 
isotropic. 

k 8nG 



.R /rV k 



d(pR3) = -pd(R3) . (2.3) 

Here R is the cosmic scale factor, p and p denote the pressure and the density, respectively, and 
G is Newton's constant. Eq. (12.11 ) is called the Friedmann equation, Eq. (I2.3D is the 1st law of 
thermodynamics. The parameter k can be chosen as ±1 or to describe spaces with constant 
positive or negative curvature, or flat geometry, respectively. Eq. (|2.1|) and Eq. (|2.2[) can be 
subtracted to yield an equation for the acceleration of the scale factor 

K 47tG , ^ , 

^ = —{p + 3p). (2.4) 

The Hubble parameter H{t) determines the expansion of the Universe. It is defined as H = ^. 
The present day value H(0) = Hq is called the Hubble constant. With this definition the critical 

density pc — the density, where the geometry of the Universe is flat — follows as pc = . To 
solve Eqs. (|2.1|) - (|2.3|) we need an additional ingredient: an equation of state, that describes the 
connection between density and pressure of the matter content of the Universe (i.e. radiation, 
baryonic matter or dark energy). At the beginning, the Universe was dominated by radiation, 
after recombination the photons decoupled and the Universe was matter dominated. Today the 
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2.1. The Cowsik-McClelland-bound 



dark energy contributes most of the density of the Universe. The equations of state are 

1 

p = -p for radiation, (2.5) 

p = for matter, (2.6) 
p = —p for dark energy. (2.7) 
The Eqs ( |2.5H2.7| ) can be summarized to 

1 

p = wp, with w = {y 0, -1}, (2.8) 

for radiation, matter, and dark energy, respectively. The dark energy is connected to the cosmo- 
logical constant in the Einstein field equation. 

2.1.2. Basic Thermodynamics 

The particle density, the energy density and the pressure of a particle species in the Universe 
are given by 



(27r)3 



(27r)3 



d3p/(p), (2.9) 
J d3p£(p)/(p), (2.10) 

/ d'p^/(p). (2.11) 



^ (271)3 

where the phase space distribution (or occupancy) is given by 

/(P) = eiE-llr^, - (2-12) 

The + sign holds for fermions, the — for bosons, and }i is the chemical potential of the particles 
species. The energy of a relativistic particle is given by E(p) = \/p^ + ni^- The entropy S 
follows from 

TdS = dpV + pdy = d[{p + p) V] - Vdp . (2.13) 
2.1.3. Particles in the Universe 

I consider the behaviour of a class of particles ^i, i = 1 . . . n (f. e. sparticles in the MSSM) in 
the thermal bath of the early Universe. Griest and Seckel discuss in [15] the mechanisms that 
are important in order to determine the number density of these new particles. They assume 
that the xpi have a multiplicatively conserved quantum number which distinguish them from 
Standard Model (SM) particles. In the MSSM, R parity [16] is such a quantum number. The 
subsequent reactions appear: 

ipiXpj ^ XX', (2.14a) 
xpiX ^ xpjX', (2.14b) 
rpj ^ tpiXX'. (2.14c) 
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2. Cosmological bounds on neutralino masses 



where X, X' denote SM particles. Examples in the MSSM for these three reaction types are: 

(2.15a) 
(2.15b) 
(2.15c) 



xVe- 



respectively. One of these particles is stabile due to the conserved quantum number. In the 
MSSM with conserved R-parity, it is the x\- Griest and Seckel classify the reaction types, see 
Eq. (|2.14|) , further. If the lightest i/?, = \p\ is nearly mass degenerate to the next to lightest particle 
Tp2, then the number density of t/^i is also determined by annihilations of ^2 which decays later 
into . This is called coannihilation. The masses of annihilation products can be heavier than 
the masses of the ingoing particles, if the energy of the ingoing particles is large enough. Griest 
and Seckel call this "forbidden" channels. If annihilation occurs at a pole in the cross section it 
is called annihilation near a pole or resonant annihilation. 

For the further discussion, I exclude coannihilation and resonant annihilation for simplicity. 

The time evolution of a particle ^ with total cross section a is described by the Boltzmann 
equation: 



dn 



df 



3Hn^ + {a\v\)\n 



^Eq] = 0, 



(2.16) 



with the the particle velocity v. The second term accounts for the dilution of the species due to 
the expansion of the Universe, the third term for the decrease by annihilation into other particles 
or coannihilation with other particles. If we define 



y 



S_ 
V 

s 

m 
Y 



P + P 
T 



{V : volume). 



(m: particle mass). 



H{m) 



m 



1.67§-y^^^ (mpi : Planck mass). 



bosons 



t = 030ig: 



fermions 



(2.17a) 
(2.17b) 
(2.17c) 

(2.17d) 

{Tj : temperature of particle species /), 

(2.17e) 
(2.17f) 



then Eq. (|2.16l l can be cast into 



dY xs , , ,> / T T \ 



or 



X dY 

YEq dx 



Li 

H 



Y 



Yi 



Eq 



nEq{0-A\ 



(2.18) 



(2.19) 
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2.1. The Cowsik-McClelland-bound 



g^, is the number of massless degrees of freedom at T,, where the particle temperature T, ac- 
counts for the possibility that it is different from the photon temperature T. The thermal aver- 
aged cross section (cr|z^|) is defined as 

{^\v\) = j-^ [n7£sk\M\\2nrS^'Hp, + p2-p3-p (2.20) 

If y = Y-Eq, then Y does not change with time, so it is constant as expected, c. f. Eq. (|2.20|) . If 
Tj^/H < 1, then the relative change of is small and the annihilation processes stop, which 
means that the number of that particle species remains constant within a comoving volume. 



2.1.4. Application to Massless Neutralinos 

In the MSSM with R-parity conservation, the lightest neutralino is stabile and can be the lightest 
supersymmetric particle. Therefore it is a dark matter candidate. I discuss the case when the 
neutralino is (nearly) massless, < 0{1 eV). The Z width allows a higgsino contribution of 

about Y^Nf3 + N^^ < 0.5 ^ (0.08)^/^, see Choudhury et al. [10]. Mi, the bino-mass, is normally 
chosen smaller than M2 and ji, and so the lightest neutralino is almost 100% bino. For simplicity, 
I assume that it is purely bino. The neutralino freezes out at Xj = m/Tj <^ 3, and at freeze out 
it is still relativistic. From that it follows Y{t — > 00) = Y-Eqixj). 



Y 



So ~2n^^^^'- ' 



(2.21) 



where ngq and s are given by Eq. (|2.9|) and (|2.17al) , sq is the present entropy density, and Z, 
denotes the Riemannian Zeta function. It is assumed that the entropy per comoving volume 
is conserved, geff counts the degrees of freedom of the neutralino field multiplied with 3/4 to 
correct for the fermionic nature of the field, g-^o, counts the number of relativistic fields at freeze 
out, whereby fermionic fields are corrected with 7/8: 



S'eff 



bosons ^ ^ ferminns ^ / 



fermions 



g, ip = boson 
|g, rp = fermion 



The neutralino density is obtained by 



Px 



m^SoY{t = 00) = 



45 ^(3) geff 

2 g^siry 



P;, _ 43 a3) 87rG geff ^3„, 



(2.22) 



(2.23) 



(2.24) 
(2.25) 



In Eq. (|2.25|) I relate the relic density Cih^ to the photon temperature by using sq = T^ and 
to the critical density. The constraint on the density is chosen such that the lightest neutralino 
does not disturb structure formation, so they cannot form the dominant component of the dark 
matter. 
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2. Cosmological bounds on neutralino masses 

Light neutralinos decouple at T 0(1 — 10 MeV). This temperature is somewhat higher 
than the temperature, when the neutrinos decouple. This is due to the selectron mass which can 
be larger than the Z mass, leading to smaller cross sections. But the temperature is below the 
muon mass so that it is not necessary to know the exact value. Nevertheless we have 2 bosonic 
and 12 fermionic relativistic degrees of freedom (one Dirac electron, three left handed neutrino 
species, one photon, one light Majorana neutralino) leading to g^^s = 12.5 and geff = l-^- H I 
demand (value of dyh^ taken from WMAP [13]) 

n^h^ < Clyh^ = 0.0067, (2.26) 

then it follows 

< 0.7/h^ eV . (2.27) 

This idea is due to Gershtein and Zel'dovich [17] and Cowsik and McClelland [18] to derive 
neutrino mass bounds. 



2.2. The Lee - Weinberg bound 

In this section, I discuss mass bounds for heavy nonrelativistic neutralinos with > 0(10 GeV) . 
I use the same method as proposed by various authors independently in [19-22] to constrain 
heavy neutrinos. This bound is now referred as Lee - Weinberg bound. 

This case is not as easy, the thermal averaged cross section and the freeze out temperature 
have to be calculated to yield an approximate solution of the Boltzmann equation. 

For simplicity, I consider only the neutralino annihilation into leptons 

X°ifi -^f^, ^ = e, ^, T, Ve, t/T. (2.28) 

The T is considered as massless0, all sleptons have common mass Mj (not to be confused with 
the common scalar mass parameter Mq), so the cross sections are related by 

<X\X\ ^ rj.t) = 16c7(;e?X? ^ = ^^'^{fiXl ^ ViVi)- (2.29) 

The thermal averaged cross section Eq. (|2.20|) can be calculated using the techniques described 
in [24]. This leads to a parametrisation of the form ~ ctqx^". In the case of a bino the 

thermal averaged cross section reads 

ml 

Hfifi - mv\) ^ = 5471—^^ x-\ (2.30) 

with X defined in Eq. I|2.17c|) . The Boltzmann Equation can be cast into the form 

dY _ f x{(7\v\)s\ 2 



-) ^x-"-2(y2(x)-y|q(x)) . (2.31) 



dx \ H{m) 

Let the difference A(x) denote the deviation Y^x) — yEq(^) of the particle density of the bino 
from equilibrium density Y^(^{x) = 0.1^5{g/ g^.s)^'^^'^^^^ ■ Shortly after the Big Bang, the de- 
viation and its derivative are small. Therefore, a good approximation is setting |^A(x)| = 
I A'(x) I 0, and one gets: 

j'l+Zy' (x) 

A(x) ^ — ^ , 1< x<t:xf. (2.32) 



^For a lower mass bound of about w 15 GeV, this is a good approximation, but not for neutralino masses of the 
order 0(1 GeV) [23]. 
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2.2. The Lee - Weinberg bound 



Later after decoupling, the neutralinos are no longer in thermal equilibrium, and the terms 
involving yEq(^) can be neglected. So one gets the following differential equation: 



A'(x) 



x{(t\v\}s 
H(m) 



.-'i-2a2 



.T = l 



(2.33) 



To solve Eq. ( I2.33D , we have to integrate from x = to x = oo. Recall, that we transformed 
the time dependence of the Boltzmann equation into an x-dependence by the transformations 
Eq. (l2.17cH2.17fD . The solutions for Eqs jzfel) and (l233l) are 



1 



A 



H{m) } 

n+1 



yn+2 



1 < X < x^, 
X| < X . 



(2.34) 



Eq. (|2.34|) requires the knowledge of the freeze out temperature Tf or, equivalently, Xf = m/Tf. 
The decoupling temperature is the temperature, when the deviation A has grown to order 
yEq(x). One sets A(x|) = cyEq(x), c = 0{1), and solves Eq. (|2.32|l for Xf, yielding 



]n[0.U5{g/gV^)in 



n + ^]ln 



y(x = oo) = A(x = 00 ) 



{cr\v\)s/H{m))^^^] 



ln[0.U5{g/gy^){n + l){x{a\v\)s/H{m)) 
3.79 (n + l)x"+^ 



x=l 



SoA(x = oo) 



1.13 X 104(n + l)x"+^ 



{g*s/gV^)mnm^o-o 
1.07X 10'^(n + l)x^+i 



cm 



-3 



{g*s/gV^)mncro 



GeV 



-1 



(2.35) 
(2.36) 
(2.37) 
(2.38) 



The choice c(c + 2) = n + 1 [14] has been implemented and yields the best fit to the relic density. 
There is a remarkable feature of Eq. ( I2.38D : The lower the cross section the larger the relic density. 
This can be understood: The particle density distribution is a Boltzmann distribution. If the 
cross section is large the particles stay longer in thermal equilibrium, and the particle density 
decreases stronger with falling temperature. 



In Fig. 2.1(a) I show contours of constant relic density in the M - m^-plane. The lower right 
hand triangle of the figure is excluded since the sleptons are lighter than the neutralino. In 
Fig. |2.1(b)[ I show the contours limiting the Odm/i^ ± Sctq = 0.113 ± 3 x 0.008 area ( [25]), 
afi denotes the absolute error on Qdm^^- The horizontal line indicates the approximate lower 
bound on the slepton masses of about 80 GeV. If we demand that the neutralinos constitute the 
whole dark matter and that the sleptons are heavier than 80 GeV, we find a lower mass bound 
of the neutralinos of about 13 GeV. The masses of the slepton cannot exceed ~ 400 GeV. If the 
next to lightest supersymmetric particle is heavier than 400 GeV, the neutralino mass bounds 



are 



13 GeV <m^< 400 GeV 



(2.39) 
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2. Cosmological bounds on neutralino masses 
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JMxf [GeV] [GeV] 

(a) Contour lines of equal relic density of a bino type (b) CIdmH^ ± 3a area of the relic density 
Isp in the - M^-plane 



Figure 2.1.: relic density of a bino type Isp 



The result shows the advantage of estimating the neutralino mass from the dark matter den- 
sity form an approximation rather than doing the full calculation: There are only two (or three) 
parameters {m^, Mj, or Mq), which can be plotted in a two dimensional figure. 

I summarize the assumptions which enter the above mass bounds (|2.39P : 

• The neutralino is a nonrelativistic bino. 

• The annihilation products are charged leptons, which are considered as massless. 

• Coannihilation and resonant annihilation is unimportant. 

• R-parity (Pg - hexality [26]) is conserved. 

2.3. Numerical solution of the full Boltzmann equation 

In the previous section, I derived from an approximate solution of the Boltzmann equation an 
upper and lower bound on the neutralino mass and - with caution - for the slepton mass. Now I 
compare these results with the exact solution. For this purpose I use the program micrOMEGAs 
[27]. 

The estimation does not take into account coannihilation and resonant annihilation. Near the 
threshold where the neutralino is almost mass degenerate with the sleptons there is coannihi- 
lation. And even for a small Higgsino component, there is large resonant annihilation if the 
neutralino mass is half of the Z'^-mass or half of the /j^-mass. 

In Fig. l2.2[ I show contour lines of the relic density for the following scenario: M2 = 200 GeV, 
}i = 300 GeV, M3 = 800 GeV, tan|6 = 10, Mu, = 450 GeV, Aj = fitan^, = 1000 GeV. The 
masses of the sleptons and of the lightest neutralino are varied. The masses of the particles other 
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2.3. Numerical solution of the full Boltzmann equation 




50 100 150 200 50 100 150 200 250 300 

[GeV] [GeV] 

(a) Qom''^ i 3(7 area of the relic density. Numerical so- (b) Fig. |2.2(a")] overlayed with Fig. |2.1(b")] to compare ap- 
lution of the full Boltzmann equation for a neutralino Isp proximate and exact solution, 
with input data: M2 = 193 GeV, }i = 350 GeV, M3 = 
800 GeV, tan^ = 10, Mh, = 450 GeV, Ar = fitanfi, 
Mq = lOOOr-'" 



10 F 




(c) Neutralino density as a function of its mass for M2 = 193 GeV, f.i = 350 GeV, 
tan ji = 10, common slepton mass = 150 GeV, common squark mass = 
1000 GeV, M3 = 800 GeV, Mh^ = 450 GeV. 



Figure 2.2.: Comparison of approximate and exact calculation of the relic density for a neu- 
tralino Isp. 
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2. Cosmological bounds on neutralino masses 



than sleptons are kept constant, so one can directly see the influence of the particle masses on 
the relic density. Electron and muon are considered as massless, and the choice of A-^ leads to 
equal slepton masses. The corner at the bottom right is excluded since the neutralino is heavier 
than the sleptons. Contrary to the Lee-Weinberg-approximation, the exact solution includes 



coannihilation near the line m^o = m|. Fig. 2.2(a) shows also the influence of a (small) Higgsino 



component, leading to resonant annihilation due to and h bosons. The resonance increases 
the cross section dramatically. This allows for larger slepton masses. The two resonances appear 



in Fig. 2.2(a) as two valleys in the ni^-Mj plane at = mz/2 and = nifj/l. From this 
I conclude that in realistic models no bound can be set on the slepton mass by relic density 
calculations. As lower bound on the neutralino mass I get 

10-15GeV<m^. (2.40) 

This agrees with the lower bound from the approximation. The upper bound is given by the 



mass of the next to lightest supersymmetric particle (nlsp). For comparison Fig. 2.2(b) shows 
the approximate and the exact solution overlayed in one plot. Apart from the valleys both plots 
agree quite well. 

If non-relativistic neutralinos constitute the whole dark matter, they cannot be completely an- 
nihilated due to resonant annihilation. This means that the mass of the neutralino is sufficiently 
far away from the relations ni^ = mz/2 or = nih/l. 



Fig. |2.2(c)| shows for one parameter set (M2 = 193 GeV, |^ = 350 GeV, tanj6 = 10, M|, = 
150 GeV as common slepton mass, Mq = 1000 GeV as common squark mass, M3 = 800 GeV, 
and = 450 GeV) the relic density Clh^ as a function of the neutralino mass. Mi has been 
increased from Mi = 1.3 GeV to 130 GeV to vary the neutralino mass. The two spikes at the end 
of the curve stem from resonant annihilation due to the and the h resonance. The nlsp has 
a mass of about 135 GeV. The qualitative shape of the curve is similar to the curve published 
in [14]. 

The horizontal red dashed lines are lines with Q/z^ = Odm'?'^ i Sc^n with Odm'?'^ = 0.113, 
(To = 0.008. The black curve crosses the allowed ribbon twice: at very light neutralinos with 
mass 0(10^^ GeV) and at massive neutralinos with mass 0(10 GeV). In the first case the parti- 
cles which constitute the dark matter cannot only be neutralinos because too many relativistic 
particles disturb structure formation in the early Universe. To avoid this constraint the neutrali- 
nos are only allowed to contribute as much as the neutrinos. This lowers the neutralino mass 
bound a little bit. The bound for relativistic neutralinos agrees very well with the predictions of 
the Cowsik-McClelland-bound. 

The exact value of the lower mass bound in the nonrelativistic case depends on the param- 
eters of the model (slepton and squark masses, mass difference to the nlsp, resonant annihila- 
tion). The upper bound is rather trivial, it is the next to lightest supersymmetric particle. Such 
searches need a lot of CPU time and have recently been done by Hooper and Plehn [28], Bottino 
and al. [23] and Belanger et al. [29]. Hooper and Plehn found a lower bound of about 18 GeV for 
a nonrelativistic neutralino, Bottino et al. found a lower bound of about 6 GeV, and Belanger et 
al. found a lower neutralino mass bound of about 6 GeV in models with a light pseudoscalar 
Higgs A with mass Ma < 200 GeV. 
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3. XrXi'P^^^^^^^^^ 



In this chapter, I derive mass bounds on the selectron mass from upper limits on the cross section 
a{e^e^ X1X2) measured by the OPAL collaboration at LEP [30], if the lightest neutralino ,^5 is 
assumed as massless. These bounds on the cross section translate into bounds on the selectron 
mass. I assume equal right and left handed selectron masses. 

The Delphi [12] and the Opal collaboration [30] have searched for SUSY particles. For neu- 
tralino pair production 

e^e- ^ xUi (3.1) 



they present upper bounds on the cross sections in the m^o-m^o plane. Their analysis assumes 

that the hadronic channels X2 ~^ ^*Xv ^* ~^ have a BR(Z* qq) = 100%. The selectron 
mass is assumed to be 500 GeV. So the two body decays into selectrons is not possible. The 
production of X1X2 ^^^^ collision occurs either by s channel exchange of a Z boson or via t and 
u channel selectron exchange. For massless neutralinos, the Xi is nearly pure bino (Nn > 0.98), 
so it couples preferably to cr, the X2 is mostly wino and couples stronger to Si. Due to the large 
selectron mass the t and u channel contributions (Te to the cross section are suppressed, so the 
dominant contribution dz comes from the s channel. The interference between Z and selectron 
exchange Cze is positive. If I denote the total cross section as az = (Tz + as + (Tzs then the e^/^ 
contribution becomes larger, if the selectron is lighter, 200 GeV < nii < 500 GeV. This ensures 
that the experimental limits on the cross section are also applicable for selectrons with mass 
< 500 GeV. Therefore, the reported bounds on the cross section are absolute upper bounds. 



In Fig. 3.1(a) I show contour lines for the cross section (7(e+e^ — > X1X2) with nii = 200 GeV 
and tan /5 = 10 in the 11-M2 plane for Mi chosen such that x\ is massless. The cross section 
reaches values up to 200 fb. There is a large parameter region where the cross section exceeds 
50 fb. From Fig. 3.1(d) taken from [30], one reads off that for a massless x\ the maximally al- 



lowed cross section is about 50 fb at a/s = 208 GeV (At m^o = 115 GeV, 125 GeV, 135 - 145 GeV, 
there are dark grey spots, indicating that the allowed cross section is 100 fb. They are most 
likely due to fluctuations in the data, I ignore them for simplicity). Within the = 104 GeV 
contour line and the 50 fb contour line the cross section is larger than 50 fb and so this part of 
the parameter space is ruled out (note that X2 ^^id x^ are nearly mass degenerate). 



In Fig. 3.1(b) I show contour lines of the minimal selectron mass so that the limits from 



Fig. 3.1(d)[ (7(e+e — > x\X2) < 50 fb, are fulfilled. The upper black line indicates the kine- 



matical limit. Below the lower black line, the Xi is lighter than 104 GeV, which is experimen- 
tally excluded [11]. Along the blue contour the X2 the selectrons have equal masses at 
about 175 GeV. Above the blue line the selectrons are lighter than xH the two body decay 
X2 ~^ ^R/L^ is allowed. For m^o > 175 GeV no part of the parameter space can be excluded. In 
Fig. 3.1(c)[ I show contour lines for the mass of X2 in the fi-Mi plane. 



For fi, M2 < 200 GeV the OPAL bound is only fulfilled if the selectrons are heavier than 
~ 350 GeV. For fi = 352 GeV, M2 = 193 GeV as in the SPSla scenario, the right handed selectron 
must be heavier than 180 GeV. 
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3. Xi-Ai'P^^^^'^^^'^^ 

Conclusion: The experimental limits on the cross section a{e^e^ X1X2) OPAL set severe 
bounds on the selectron mass if m^o =0 GeV and m^o < 175 GeV. 

Ai A2 




Observed limit at 208 GeV, Hadronic Channel, z" BR 




(c) Contour lines for the X2 mass in GeV in the pi - (d) Observed 95% confidence level upper bound on the cross 
M2 plane section <7(e+e^ based on decays to hadronic final 

states and assuming 100% branching ratio for decay into Z"*. 

This figure is taken from [30]. 

Figure 3.1.: Deriving a lower mass bound on the selectron mass if Xi is massless. 
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4. Radiative Neutralino Production 



4.1. Introduction 

Supersymmetry (SUSY) is an attractive concept for theories beyond the Standard Model (SM) of 
particle physics. SUSY models like the Minimal Supersymmetric Standard Model (MSSM) [1, 
5, 31] predict SUSY partners of the SM particles with masses of the order of a few hundred 
GeV. Their discovery is one of the main goals of present and future colliders in the TeV range. 
In particular, the international e+e^ linear collider (ILC) will be an excellent tool to determine 
the parameters of the SUSY model with high precision [32-36]. Such a machine provides high 
luminosity L = 500 fb^^, a center-of-mass energy of a/s = 500 GeV in the first stage, and a 
polarised electron beam with the option of a polarised positron beam [37]. 

The neutralinos are the fermionic SUSY partners of the neutral gauge and CP-even Higgs 
bosons. Since they are among the lightest particles in many SUSY models, they are expected to 
be among the first states to be observed. At the ILC, they can be directly produced in pairs 

e^^e- ^t^f], (4.1) 

which proceeds via Z boson and selectron exchange [38,39]. At tree level, the neutralino sector 
depends only on the four parameters Mi, M2, }i, and tanj6, which are the L[(l)y and SU{2)l 
gaugino masses, the higgsino mass parameter, and the ratio of the vacuum expectation values 
of the two Higgs fields, respectively. These parameters can be determined by measuring the 
neutralino production cross sections and decay distributions [35, 40-43]. In the MSSM with 
R-parity (or proton hexality Pg, [26]) conservation, the lightest neutralino Xi is typically the 
lightest SUSY particle (LSP) and as such is stable and a good dark matter candidate [44,45]. In 
collider experiments the LSP escapes detection such that the direct production of the lightest 
neutralino pair 

e^e- ^ xlxl (4.2) 

is invisible. Their pair production can only be observed indirectly via radiative production 

e^e- ^ xUh, (4.3) 

where the photon is radiated off the incoming beams or off the exchanged selectrons. Although 
this higher order process is suppressed by the square of the additional photon-electron cou- 
pling, it might be the lightest state of SUSY particles to be observed at colliders. The signal is a 
single high energetic photon and missing energy, carried by the neutralinos. 

As a unique process to search for, the first SUSY signatures at e+e^ colliders, the radiative 
production of neutralinos has been intensively studied in the literature [46-62] |^ Early investi- 
gations focus on LEP energies and discuss special neutralino mixing scenarios only, in particular 
the pure photino case [46-53]. More recent studies assume general neutralino mixing [54-62] 
and some of them underline the importance of longitudinal [54-57] and even transverse beam 



^In addition I found two references [63,64], which are however almost identical in wording and layout to Ref. [58]. 
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4. Radiative Neutralino Production 

polarisations [54,57]. The transition amplitudes are given in a generic factorised form [55], 
which allows the inclusion of anomalous WW7 couplings. Cross sections are calculated with 
the program CompHEP [56], or in the helicity formalism [57]. Some of the studies [58-62] however 
do not include longitudinal beam polarisations, which might be essential for measuring radia- 
tive neutralino production at the ILC. Special scenarios are considered, where besides the sneu- 
trinos also the heavier neutralinos [59-61], and even charginos [65-67] decay invisibly or almost 
invisibly. However, a part of such unconventional signatures are by now ruled out by LEP2 
data [59,66,68]. For the ILC, such "effective" LSP scenarios have been analysed [60], and strate- 
gies for detecting invisible decays of neutralinos and charginos have been proposed [65, 67]. 
Moreover, the radiative production of neutralinos can serve as a direct test to see, whether 
neutralinos are dark matter candidates. See for example Ref. [69], which presents a model in- 
dependent calculation for the cross section of radiatively produced dark matter candidates at 
high-energy colliders, including polarised beams for the ILC. 

The signature "photon plus missing energy" has been studied intensively by the LEP collab- 
orations ALEPH [70], DELPHI [71], L3 [72], and OPAL [68,73]. In the SM, 

e^e^ vvj (4.4) 

is the leading process with this signature. Since the cross section depends on the number Ny 
of light neutrino generations [74], it has been used to measure Np consistent with three. In 
addition, the LEP collaborations have tested physics beyond the SM, like non-standard neutrino 
interactions, extra dimensions, and SUSY particle productions. However, no deviations from 
SM predictions have been found, and only bounds on SUSY particle masses have been set, e.g. 
on the gravitino mass [70-73]. This process is also important in determining collider bounds on 
a very light neutralino [75]. For a combined short review, see for example Ref. [76]. 

Although there are so many theoretical studies on radiative neutralino production in the lit- 
erature, a thorough analysis of this process is still missing in the light of the ILC with a high 
center-of-mass energy, high luminosity, and longitudinally polarised beams. As noted above, 
most of the existing analyses discuss SUSY scenarios with parameters which are ruled out by 
LEP2 already, or without taking beam polarisations into account. In particular, the question of 
the role of the positron beam polarisation has to be addressed. If both beams are polarised, the 
discovery potential of the ILC might be significantly extended, especially if other SUSY states 
like heavier neutralino, chargino or even slepton pairs are too heavy to be produced at the 
first stage of the ILC at i/s = 500 GeV. Moreover, the SM background photons from radiative 
neutrino production e+e^ — > vv'y have to be included in an analysis with beam polarisations. 
Proper beam polarisations could enhance the signal photons and reduce those from the SM 
background at the same time, which enhances the statistics considerably. In this respect also the 
MSSM background photons from radiative sneutrino production 

e^e^ vv*j (4.5) 

have to be discussed, if sneutrino production is kinematically accessible and if the sneutrino 
decay is invisible. 

Finally, the studies which analyse beam polarisations do not give explicit formulas for the 
squared matrix elements, but only for the transition amplitudes [54,55,57]. Other authors admit 
sign errors [61] in some interfering amplitudes in precedent works [60], however do not provide 
the corrected formulas for radiative neutrino and sneutrino production. Additionally, I found 
inconsistencies and sign errors in the Z exchange terms in some works [54, 57], which yield 
wrong results for scenarios with dominating Z exchange. Thus I will give the complete tree- 
level amplitudes and the squared matrix elements including longitudinal beam polarisations. 
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Figure 4.1.: Diagrams for radiative neutralino production e^e^ XiXil [77]. For the calcula- 
tion in Appendix |Al the first number of the diagrams labels f-channel, the second 
one M-channel exchange of selectrons, where the neutralinos are crossed. 



such that the formulas can be used for further studies on radiative production of neutralinos, 
neutrinos and sneutrinos. 

In Sec. I4.2[ 1 discuss my signal process, radiative neutralino pair production, as well as the 
major SM and MSSM background processes. In Sec. l4.3[ 1 define cuts on the photon angle and 
energy, and define a statistical significance for measuring an excess of photons from radiative 
neutralino production over the backgrounds. I analyse numerically the dependence of cross 
sections and significances on the electron and positron beam polarisations, on the parameters of 
the neutralino sector, and on the selectron masses. I summarise and conclude in Sec. 14.51 In the 
Appendix, I define neutralino mixing and couplings, and give the tree-level amplitudes as well 
as the squared matrix elements with longitudinal beam polarisations for radiative production 
of neutralinos, neutrinos and sneutrinos. In addition, I give details on the parametrisation of 
the phase space. 



4.2. Radiative Neutralino Production and Backgrounds 
4.2.1. Signal Process 

Within the MSSM, radiative neutralino production [46-62] 

e+ +e- ^ Xi + xi + j (4.6) 

proceeds at tree-level via t- and w-channel exchange of right and left selectrons Sr^i, as well as 
Z boson exchange in the s-channel. The photon is radiated off the incoming beams or the ex- 
changed selectrons; see the contributing diagrams in Fig. 14.11 I give the relevant Feynman rules 
for general neutralino mixing, the tree-level amplitudes, and the complete analytical formulas 
for the amplitude squared, including longitudinal electron and positron beam polarisations, in 
Appendix [a] I also summarise the details of the neutralino mixing matrix there. For the calcu- 
lation of cross sections and distributions I use cuts, as defined in Eq. (|4.13|) . An example of the 
photon energy distribution and the \/s dependence of the cross section is shown in Fig.[ 
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4.2.2. Neutrino Background 

Radiative neutrino production [60, 74, 78-80] 

+ vg + vi + j , 1 = e,}i,T (4.7) 

is a major SM background. Electron neutrinos Ve are produced via f-channel W boson exchange, 
and V£',;i,T via s-channel Z boson exchange. I show the corresponding diagrams in Appendix [HI 
where I also give the tree-level amplitudes and matrix elements squared including longitudinal 
beam polarisations. 

4.2.3. MSSM Backgrounds 

Next 1 consider radiative sneutrino production [60, 81, 82] 

+ ^ V£ + i7;+7, l = e,}i,T. (4.8) 

I present the tree-level Feynman graphs as well as the amplitudes and amplitudes squared, in- 
cluding beam polarisations, in Appendix|C] The process has f-channel contributions via virtual 
charginos for i/gV* -production, as well as s-channel contributions from Z boson exchange for 
i?e,/(,Tf'e ;^,T"Production, see Fig. IC.ll Radiative sneutrino production, Eq. (|4.8|) , can be a major 
MSSM background to neutralino production, Eq. (|4.6|), if the sneutrinos decay mainly invisi- 
bly, e.g., via V —>■ Xi'^- This leads to so called "virtual LSP" scenarios [59-61]. However, if 
kinematically allowed, other visible decay channels like v — > X^i^ reduce the background 
rate from radiative sneutrino production. For example in the SPS la scenario [83, 84], 1 have 
BR{ve Xi^e) = 85%, see TableHU 

In principle, also neutralino production e+e^ X1X.2 followed by the subsequent radiative 
neutralino decay [85] X2 ~^ Xil is a potential background. However, significant branching 
ratios BR(;^2 ~^ Xil) > 10% ^re only obtained for small values of tan/3 < 5 and/or Mi ~ 
M2 [62,86,87]. Thus I neglect this background in the following. For details see Refs. [86-88]. 

4.3. Numerical Results 

I present numerical results for the tree-level cross section for radiative neutralino production, 
Eq. (|4.6|) , and the background from radiative neutrino and sneutrino production, Eqs. (|4.7|) and 
(|4.8|) , respectively. I define the cuts on the photon energy and angle, and define the statistical 
significance. 1 study the dependence of the cross sections and the significance on the beam po- 
larisations Pf,- and Pg+, the supersymmetric parameters fi and M2, and on the selectron masses. 
In order to reduce the number of parameters, I assume the SUSY GUT relation 

Ml = ^ tan^ 0,„M2 . (4.9) 

Therefore the mass of the lightest neutralino is m^o > 50GeV [89]. I also use the approximate 
renormalisation group equations (RGE) for the slepton masses [90-92], 

m?^ = ml + 0.23Ml-mlcos2^sin^e,u, (4.10) 
m?^ = m^ + O.79M| + m|cos2,6(^-^ + sin2 0j„), (4.11) 

ml = m§ + 0.79M^ + ^m|cos2(6, (4.12) 

with niQ the common scalar mass parameter. Since in my scenarios the dependence on tan j6 is 
rather mild, I fix tan j6 = 10. 
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Table 4.1.: Parameters and masses for SPS la scenario [83,84]. 



tan j6 = 10 


= 352 GeV 


M2 = 193 GeV 


mo = 100 GeV 


m^o = 94 GeV 


m^± = 178 GeV 


mg^ = 143 GeV 


niy^ = 188 GeV 


m n = 178 GeV 

A7 


m^± = 376 GeV 

A2 


niei^ = 204 GeV 


BR(i7, ^ = 85% 



4.3.1. Cuts on Photon Angle and Energy 

To regularise the infrared and collinear divergencies of the tree-level cross sections, I apply cuts 
on the photon scattering angle and on the photon energy Ey 

-0.99 < cos0^ < 0.99, 0.02 < x < 1 x = 7—^, (4.13) 

^beam ^beam 

with the beam energy Ebeam — \/s/2. The cut on the scattering angle corresponds to 9^ S 
[8°, 172°], and reduces much of the background from radiative Bhabha scattering, e+e^ — ^ 
e+e^7, where both leptons escape close to the beam pipe [70, 71]. The lower cut on the pho- 
ton energy is E-y = 5 GeV for -^s = 500 GeV. The upper cut on the photon energy x""^" = 
1 — m^o / ^beam kinematical limit of radiative neutralino production. At i/s = 500 GeV 

and for m^o > 70 GeV, this cut reduces much of the on-shell Z boson contribution to radia- 
tive neutrino production, see Refs. [56,59,82,93] and Fig. |4.2[ a). I assume that the neutralino 
mass m^o is known from LHC or ILC measurements [35]. If m^o is unknown, a fixed cut, e.g., 
E^ax ^ at ^ = 500 GeV, could be used instead [93]. 

4.3.2. Theoretical Significance 

In order to quantify whether an excess of signal photons from neutralino production, Ng = aC, 
for a given integrated luminosity C, can be measured over the SM background photons, Nb = 
a^C, from radiative neutrino production, I define the theoretical significance S and the signal to 
background ratio r (or reliability) 



'^Signal 
•^Background 



(4.14) 
(4.15) 



A theoretical significance of, e.g., S = 1 implies that the signal can be measured at the statistical 
68% confidence level. Also the the signal to background ratio Ns/Nb should be considered 
to judge the reliability of the analysis. For example, if the background cross section is known 
experimentally to 1% accuracy, I should have Ns/Nb > 1/100. 

I will not include additional cuts on the missing mass or on the transverse momentum distri- 
butions of the photons [56,93]. Detailed Monte Carlo analyses, including detector simulations 
and particle identification and reconstruction efficiencies, would be required to predict the sig- 
nificance more accurately, which is however beyond the scope of the present work. Also the 
effect of beamstrahlung should be included in such an experimental analysis [93-95]. Beam- 
strahlung distorts the peak of the beam energy spectrum to lower values of Ebeam = \/s/2, and 
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Figure 4.2.: (a) Photon energy distributions for a/s = 500 GeV, and (b) ^/s dependence of 
the cross sections cr for radiative neutralino production e+e^ XiXil (black, 
solid), neutrino production e+e^ vvj (violet, dashed) and sneutrino produc- 
tion e+e^ {/{/*7 (blue, dotted) for scenario SPS la [83,84], see Table |4T1 with 
{Pe-, Pe+) = (0.8, —0.6). The red dot-dashed line is in (a) the photon energy distri- 
bution for radiative neutrino production e^e^ vv^, and in (b) the cross section 
without the upper cut on the photon energy £-y, see Eq. (|4.13p . 

is more significant at colliders with high luminosity. In the processes I consider, the cross sec- 
tions for e+e^ — > and e+e^ vvj depend significantly on the beam energy only near 
threshold. In most of the parameter space we consider, for a/s = 500 GeV the cross sections are 
nearly constant, see for example Fig. I4.2r b), so I expect that the effect of beamstrahlung will be 
rather small. However, for M2, }i > 300 GeV, e+e^ XiXil is the only SUSY production pro- 
cess, which is kinematically accessible, see Fig. 14.41 In order to exactly determine the kinematic 
reach, the ILC beamstrahlung must be taken into account. 

4.3.3. Energy Distribution and ^/s Dependence 

In Fig. I4.2f a), I show the energy distributions of the photon from radiative neutralino pro- 
duction, neutrino production, and sneutrino production for scenario SPS la [83, 84], see Ta- 
ble |4]TJ with ^/s = 500 GeV, beam polarisations (Pg-,Pg+) = (0.8,-0.6), and cuts as defined 
in Eq. (|4.13|) . The energy distribution of the photon from neutrino production peaks at Ey = 
(s — m|) / (2a/s) 242 GeV due to radiative Z return, which is possible for ^/s > niz- Much of 
this photon background from radiative neutrino production can be reduced by the upper cut on 
the photon energy x"^^^ = E^'^'^/Ebeam = 1 ~ / ^beam' ^"l- SilSll/ which is the kinematical 
endpoint E™*^ w 215 GeV of the energy distribution of the photon from radiative neutralino 
production, see the solid line in Fig. l4.2f a). Note that in principle the neutralino mass could be 
determined by a measurement of this endpoint EJ^"^" = EJj"^''(m^o) 

ml, = \{s-l^sE^^-). (4.16) 

For this one would need to be able to very well separate the signal and background processes. 
This might be possible if the neutralino is heavy enough, such that the endpoint is sufficiently 
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removed from the Z'^-peak of the background distribution. 

In Fig. I4.2r b) I show the a/s dependence of the cross sections. Without the upper cut on 
the photon energy x™^", see Eq. (14.131) , the background cross section from radiative neutrino 
production e+e^ — > vv'y, see the dot-dashed line in Fig. |4.2[ b), is much larger than the corre- 
sponding cross section with the cut, see the dashed line. However with the cut, the signal cross 
section from radiative neutralino production, see the solid line, is then only about one order of 
magnitude smaller than the background. 

4.3.4. Beam Polarisation Dependence 

In Fig. I4.3f a) 1 show the beam polarisation dependence of the cross section cr{e^e^ XiXil) 
for the SPS la scenario [83,84], where radiative neutralino production proceeds mainly via right 
selectron exchange. Since the neutralino is mostly bino, the coupling to the right selectron 
is more than twice as large as to the left selectron. Thus the contributions from right selectron 
exchange to the cross section are about a factor 16 larger than the ei contributions. In addition 
the ei contributions are suppressed compared to the Sr contributions by a factor of about 2 since 
niif, < nigi^, see Eqs. (|4.10|) - (|4.11|) . The Z boson exchange is negligible. The background process, 
radiative neutrino production, mainly proceeds via W boson exchange, see the corresponding 
diagram in Fig. IB.ll Thus positive electron beam polarisation Pg- and negative positron beam 
polarisation Pg+ enhance the signal cross section and reduce the background at the same time, 
see Figs. l4.3f a) and l4.3f c), which was also observed in Refs. [55,93]. The positive electron beam 
polarisation and negative positron beam polarisation enhance exchange and suppress 
exchange, such that it becomes negligible. Opposite polarisations would lead to comparable 
contributions from both selectrons. In going from unpolarised beams {Pe-,Pe+) = (0,0) to 
polarised beams, e.g., (Pg-,Pg+) = (0.8, —0.6), the signal cross section is enhanced by a factor 
~ 3, and the background cross section is reduced by a factor 10. The signal to background 
ratio increases from Ns/Nq ~ 0.007 to Nc,/Nb ~ 0.2, such that the statistical significance S, 
shown in Fig. I4.3f b), is increased by a factor 8.5 to S ~ 77. If only the electron beam is 
polarised, {Pe-,Pe+) = (0.8,0), I still have Ns/Nb ^ 0.06 and S ^ 34, thus the option of a 
polarised positron beam at the ILC doubles the significance for radiative neutralino production, 
but is not needed or essential to observe this process at ^/s = 500 GeV and C = 500 fb^^ for the 
SPS la scenario. 

In contrast, the conclusion of Ref. [56] is, that an almost pure level of beam polarisations 
is needed at the ILC to observe this process at all. The authors have used a scenario with 
Ml = Ml, leading to a lightest neutralino, which is mostly a wino. Thus larger couplings 
to the left selectron than to the right selectron are obtained. In such a scenario, one cannot 
simultaneously enhance the signal and reduce the background. Moreover their large selectron 
masses m^^j^ = 500 GeV lead to an additional suppression of the signal, see also Sec. 14.3.6] 

Finally I note that positive electron beam polarisation and negative positron beam polarisa- 
tion also suppress the cross section of radiative sneutrino production, see Fig. l4.3f d). Since it is 
the corresponding SUSY process to radiative neutrino production, I expect a similar quantita- 
tive behaviour. 

4.3.5. }i &. M2 Dependence 

In Fig. l4.4f a), I show contour lines of the cross section C7"(e+e^ XiXij) in fb in the (f/, M2)- 
plane. For fi > 300 GeV the signal and the background cross sections are nearly independent of 
ji, and consequently also the significance, which is shown in Fig. l4.4r b). In addition, the depen- 
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Figure 4.3.: (a) Contour Imes of the cross section and (b) the significance S for e^e — > XiXil ^f 
a/s = 500 GeV and £. = 500 fb"^ for scenario SPS la [83,84], see TableHH The beam 
polarisation dependence of the cross section for radiative neutrino and sneutrino 
production are shown in (c) and (d), respectively. 
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Figure 4.4.; Contour lines (solid) of (a) the cross section a{e^e XiXil)' (b) the significance S, 
(c) the neutrino background cr^ (e+e^ — > 1/1/7), ^rid (d) the neutralino mass m^o in the 

^-M2 plane for = 500 GeV, (Pe-,Pg+) = (0.8, -0.6), C = 500 fb"\ with tan ,6 = 
10, mo = 100 GeV, and RGEs for the selectron masses, see Eqs. (I4.10D , ( 14.111) . The 
grey area is excluded by m^± < 104 GeV. The dashed line indicates the kinematical 
limit m^o + m^o = \/s, and the dot-dashed line the kinematical limit 2m ^± = \/s. 
Along the dotted line in (b) the signal to background ratio is t7/(7B = 0.01. The area 
A is kinematically forbidden by the cut on the photon energy E^, see Eq. 04.131 1. 
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dence of the neutralino mass on }i is fairly weak for > 300 GeV, as can be seen in Fig. |4.4f d). 
Also the couplings have a rather mild /^-dependence in this parameter region. 

The cross section a^{e^e^ — > 1/1/7) of the SM background process due to radiative neutrino 
production, shown in Fig. I4.4r c), can reach more than 340 fb and is considerably reduced due 
to the upper cut on the photon energy x"^^^, see Eq. (|4.13|) . Without this cut I would have erg = 
825 fb. Thus the signal can be observed with high statistical significance S, see Fig. l4.4r b). Due 
to the large integrated luminosity C = 500 fb"^ of the ILC, I have S > 25 with Ns/Nb > 1/4 for 
M2 < 350 GeV. For }i < 0,1 get similar results for the cross sections in shape and size, since the 
dependence of Nn on the sign of fi, see Eq. (|A.3|) , is weak due to the large value of tan j6 = 10. 

In Fig. 14.41 I also indicate the kinematical limits of the lightest observable associated neu- 
tralino production process, e+e^ X1X2 (dashed), and those of the lightest chargino produc- 
tion process, e+e^ XiXi (dot-dashed). In the region above these lines M2 > 300 GeV, 
heavier neutralinos and charginos are too heavy to be pair-produced at the first stage of the 
ILC with ^/s = 500 GeV. In this case radiative neutralino production e+e^ X^X^l will be 
the only channel to study the gaugino sector. Here significances of 5 < S < 25 can be ob- 
tained for 350 GeV < M2 < 450 GeV, see Fig.SiHb). Note that the production of right sleptons 
e+e^ ^R^R' ^ = ^/ srid in particular the production of the lighter staus e+e^ ^i^^i / due 
to mixing in the stau sector [96], are still open channels to study the direct production of SUSY 
particles for M2 < 500 GeV in our GUT scenario with mo = 100 GeV. 

4.3.6. Dependence on the Selectron Masses 

The cross section for radiative neutralino production cr{e^e^ XiXil) proceeds mainly via se- 
lectron gR L exchange in the t and w-channels. Besides the beam polarisations, which enhance Sr 
or ei exchange, the cross section is also very sensitive to the selectron masses. In the mSUGRA 
universal supersymmetry breaking scenario [97], the masses are parametrised by niQ and M2, 
besides tan fi, which enter the RGEs, see Eqs. ([4.10|) and l|4.11|l . I show the contour lines of the 
selectron masses Cr^l in the mo-M2 plane in Fig. I4.5f c) and I4.5f d), respectively. The selectron 
masses increase with increasing niQ and M2. 

For the polarisations (P^-, Pe+) = (0.8, —0.6), the cross section is dominated by Sr exchange, 
as discussed in Sec. 14.3.4] In Fig. l4.5f a) and l4.5f b), I show the wq and M2 dependence of the cross 
section and the significance S, Eq. (|4.14|l . With increasing mo and M2 the cross section and the 
significance decrease, due to the increasing mass of Sr. In Fig. l4.4f d). I see that for fi > 7/ 10 M2, 
the neutralino mass ra^o is practically independent of }i and rises with M2. Thus for increasing 

M2, and thereby increasing neutralino mass, the cross section cr{e^e^ XiXil) reaches the 
kinematical limit at M2 ~ 500 GeV for i/s = 500 GeV. A potential background from radiative 
sneutrino production is only relevant for M2 < 200 GeV, nio < 200 GeV. For larger values the 
production is kinematically forbidden. 

In Fig. I4.5[ I also indicate the kinematical limit of associated neutralino pair production m^o + 
m^o = y/s = 500 GeV, reached for M2 « 350 GeV. If in addition nio > 200 GeV, also selectron 

A2 

and smuon pairs cannot be produced at ^/s = 500 GeV due to mj,^ > 250 GeV. Thus, in this 
parameter range where M2 > 350 GeV and mo > 200 GeV, radiative production of neutralinos 
will be the only possible production process of SUSY particles, if I neglect stau mixing. A statis- 
tical significance of S > 1 can be obtained for selectron masses not larger than niej^ ~ 500 GeV, 
corresponding to mo < 500 GeV and M2 < 450 GeV. Thus radiative neutralino production ex- 
tends the discovery potential of the ILC in the parameter range mo G [200,500] GeV and 
M2 G [350,450] GeV. Here, the beam polarisations will be essential, see Fig. 14.61 I show con- 
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Figure 4.5.: (a) Contour lines of the cross section t7(e+e XiXil)' (b) the significance S, and 
(c), (d) the selectron masses nig^^, mg^, respectively, in the mo-M2 plane for a/s = 
500 GeV, {Pe-,Pe+) = (0.8, -0.6), £, = 500 fb"\ with }i = 500 GeV, tan^ = 10, and 
RGBs for the selectron masses, see Eqs. (|4.10|l , (|4.1H) . The dashed line indicates the 
kinematical limit m^o + m^o = ^/s. The grey area is excluded by m^± < 104 GeV, the 
area A is kinematically forbidden. 
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Figure 4.6.: Contour lines of the significance S in the niQ-M2 plane for different beam polar- 
isations {Pe-,Pe+) = (0.8,-0.6) (solid), (P,-,Pe+) = (0.8,0) (dot-dashed), and 
{P,-,P,+ ) = (0,0) (dotted), for ^ = 500 GeV, C = 500 fb"\ ^ = 500 GeV, 
tan /3 = 10, and RGBs for the selectron masses, see Eqs. (|4.10|) . (|4.1Hl . The dashed 
line indicates the kinematical limit m^o + /n^o = \/s. The grey area is excluded by 
m ± < 104 GeV, the area A is kinematically forbidden. 



tour lines of the statistical significance S for three different sets of (Pg-,Pe+). The first set has 
both beams polarised, (Pe-,Pe+) = (0.8, —0.6), the second one has only electron beam polari- 
sation, {Pe-,Pe+) = (0.8,0), and the third has zero beam polarisations {Pe-,Pe+) = (0,0). The 
beam polarisations significantly enhance the discovery potential of the ILC. At least electron 
polarisation Pg- = 0.8 is needed to extend an exploration of the mQ-M2 parameter space. 



4.3.7. Note on LEP2 

I have also calculated the unpolarised cross sections and the significances for radiative neu- 
tralino production at LEP2 energies = 200 GeV, for a luminosity of C = 100pb^\ I have 
used the cuts | cos0-y| < 0.95 and 0.2 < x < 1 — w^^o/^beam' ' I'^-ISI ). Even for rather small 

selectron masses nigj^^^ = 80 GeV, the cross sections are not larger than 100 fb. Even if I alter 
the GUT relation, Eq. (|49l) , to Mi = ri2 M2, and vary ri2 within the range 0.01 < ri2 < 0.5, I 
only obtain statistical significances of S < 0.2. These values have also been reported by other 
theoretical studies at LEP2 energies, see for example Ref. [62]. 
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If I drop the GUT relation. Mi is a free parameter. For 

^ M2m|siri(2^)sin^g^ 
^ ^M2 - m| sin(2(6) cos2 0j„ 

the neutralino is massless [9] at tree-level and is apparently experimentally allowed [75]. A 
massless neutralino should enlarge the cross section for radiative neutralino production due to 
the larger phase space, although the coupling is also modified to almost pure bino. However, 
I still find S = ©(lO^^) at most. This is in accordance with the experimental SUSY searches in 
photon events with missing energy at LEP [68, 70-73], where no evidence of SUSY particles was 
found. 



4.4. The Role of Beam polarization for Radiative Neutralino 
Production at the ILC 

4.4.1. Introduction 

Detailed measurements of the masses, decay widths, couplings, and spins of the discovered 
particles are only possible at the international linear collider (ILC) [32-35]. In the first stage of 
the ILC, the center-of-mass energy will be = 500 GeV and the luminosity, C, will be 500 fb^^ 
per year. 

In preparing for the ILC, there is an on-going debate over the extent of beam polarization to be 
included in the initial design [37,98-100]. It is clear that there will be at least 80% polarization of 
the electron beam, possibly even 90% [101]. A polarized positron beam is technically and finan- 
cially more involved. However, it is possible to achieve 30% polarization already through the 
undulator based production of the positrons [100]. In light of this discussion, it is the purpose of 
this section to reconsider the effect of various degrees of electron and positron polarization on a 
particular supersymmetric production process, namely the radiative production of the lightest 
neutralino mass eigenstate Xi 

e^ + e- ^fi+xl + r (4.18) 

I shall focus on specific regions of the supersymmetric parameter space. The signal is a single, 
highly energetic photon and missing energy, carried by the neutralinos. 

The process (|4.18|l was previously studied within the MSSM and with general neutralino 
mixing in Refs. [59-62]. The additional effect of polarized beams was considered in Refs. [55, 
56, 102]. In Ref. [102], it was shown that polarized beams significantly enhance the signal and 
simultaneously suppress the main SM photon background from radiative neutrino production, 

e^ + e'^v + v + j. (4.19) 

Moreover, it was pointed out that for certain regions of the MSSM parameter space, the pro- 
cess (|4.18|l is kinematically the only accessible SUSY production mechanism in the first stage of 
the ILC at ^/s = 500 GeV [102]. Here the heavier electroweak gauginos and the sleptons are too 
heavy to be pair produced, i.e. their masses are above 250 GeV. 

Other than the standard center-of-mass energy, ^/s = 500 GeV, at the ILC, also lower energies 
are of particular interest, namely for Higgs and top physics. Higgs strahlung, 

e+ +e- ^Z + K (4.20) 
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can be well studied at the threshold energy \/s = m;, + niz, which is a/s ~ 220 GeV, for a Higgs 
boson mass of mjj ^ 130 GeV. The CP-quantum number and the spin of the Higgs boson can be 
determined from an energy scan of the production cross section near the threshold [103]. 

From a scan at the threshold energy of top pair production, ^/s = Init ~ 350 GeV, the top 
mass mt can be determined with an error Snit < 0.1 GeV [104]. Thus the present error on the top 
mass, Snit « 3 GeV [11], and the foreseen error from LHC measurements, Smt ~ IGeV [105], 
can be reduced by one order of magnitude. Also the top width, Tf, and the strong coupling 
constant, a,/ can be precisely determined by a multi parameter fit of the cross section, top mo- 
mentum distribution, and forward-backward charge asymmetry near threshold [106]. 

In this section, I take these physics questions as a motivation to study the role of polarized 
beams in radiative neutralino production at the energies \/s = 220 GeV, 350 GeV, and 500 GeV 
at the ILC. For each beam energy, I shall focus on a specific supersymmetric parameter set 
within the context of minimal supergravity grand unification (mSUGRA) [107]. I thus consider 
three mSUGRA scenarios, which I label A, B and C, respectively, and which are listed below in 
Table |4!2| together with the resulting spectra in Table [431 1 restrict myself to mSUGRA scenarios, 
in order to reduce the number of free parameters and since I find it suffcient to illustrate my 
point. The specific scenarios are chosen such that radiative neutralino production is the only 
supersymmetric production mechanism which is kinematically accessible at the given center-of- 
mass energy. It is thus of particular interest to learn as much about supersymmetry as is possible 
through this mechanism. As I shall see, beam polarization is very helpful in this respect. 

In Sect. I4.3.1[ I define the significance, the signal to background ratio and define a first set 
of experimental cuts. In Sect. 14.4.21 I study numerically the dependence of the signal cross 
section and the SM background, the significance, and the signal to background ratio on the 
beam polarization. In particular, I compare the results for different sets of beam polarizations, 
{P,+ \P,-) - ( 0|0), (0|0.8), (-0.3|0.8), (-0.6|0.8), (0|0.9) and (-0.3|0.9). I summarize and con- 
clude in Sect. (Ill 

4.4.2. Numerical results 

I choose the three scenarios in such a way, that only the lightest neutralinos can be radiatively 
produced for each of the \/s values, respectively. The other SUSY particles, i.e. the heavier neu- 
tralinos and charginos, as well as the sleptons and squarks are too heavy to be pair produced 
at the ILC. It is thus of paramount interest to have an optimal understanding of the signa- 
ture (|4.18|) , in order to learn as much as possible about SUSY at a given ILC beam energy. Note 
that in the three scenarios (A,B,C) the squark and gluino masses are below {600, 800, 1000} GeV, 
respectively and should be observable at the LHC [35]. 

Scenario A is related to the Snowmass point SPSla [36, 83, 84] by scaling the common scalar 
mass Mo, the unified gaugino mass M1/2, and the common trilinear coupling Aq by 0.9. Thus 
the slope Mq = —Aq = 0.4 Mi /2 remains unchanged. For scenarios B and C, I also choose 
Mo = —Aq, however I change the slopes to Mq = 0.42 M1/2 in scenario B, and Mq = 0.48 M1/2 
in scenario C. For all scenarios, I fix the ratio tan /3 = 10 of the vacuum expectation values of 
the two neutral Higgs fields. In Table I4.2[ I explicitly give the relevant low energy mSUGRA 
parameters for all scenarios. These are the U{1) and SU{2) gaugino mass parameters Mi and 
M2, respectively, and the Higgsino mass parameter fi. The masses of the light neutralinos, 
charginos, and sleptons are given in Table 14.31 All parameters and masses are calculated at 
one-loop order with the computer code SPheno [108]. 
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Note that the lightest neutralino, is mostly bino in all three scenarios; 98% in scenario 
A, 99.1% in scenario B, and 99.5% in scenario C. Thus in my scenarios, radiative neutralino 
production proceeds mainly via right selectron exchange in the t and u charmel. Left selec- 
tron exchange and Z boson exchange are severely suppressed [102]. The background process 
e+e^ — > 1/1/7 mainly proceeds via W boson exchange. Thus positive electron beam polarization 
Pg- > and negative positron beam polarization Pg+ < should enhance the signal rate and 
reduce the background at the same time [55, 102]. This effect is clearly observed in Figs. l4.7[|4.8[ 
and 14.91 for all scenarios. The signal cross section and the background vary by more than one 
order of magnitude over the full polarization range. 



Table 4.2.: Definition of the mSUGRA scenarios A, B, and C. All values are given in GeV. I have fixed 
tan /3 = 10. For completeness I have included the corresponding value of \/s for each sce- 
nario. 



scenario 




Mo 


Ml/2 


Ao 


Ml 


M2 


F 


A 


220 


90 


225 


-90 


97.5 


188 


316 


B 


350 


135 


325 


-135 


143 


272 


444 


C 


500 


200 


415 


-200 


184 


349 


560 



Table 4.3.: Spectrum of the lighter SUSY particles for scenarios A, B, and C, calculated with SPheno [108]. 

All values are given in GeV. For completeness I have included the corresponding value of y/s 
for each scenario. 



scenario 






m n 

A.2 


m ± 








niy 


A 


220 


92.4 


172 


172 


124 


133 


189 


171 


B 


350 


138 


263 


263 


183 


191 


270 


258 


C 


500 


180 


344 


344 


253 


261 


356 


347 



For scenario A, I show the beam polarization dependence of the signal cross section cr{e^e^ — > 
XiXil) Fig. 4.7(a) and the dependence of the background cross section a{e'^e^ — ^ 1/1/7) iri 
Fig. 4.7(b) In both cases I have implemented the cuts of Eq. (|4.13|) . The cont lines in the Pg--Pg+ 
plane of the significance S, Eq. (|4.14|l , and the signal to background ratio r, Eq. (|4.15|l , are shown 

and those 



in Figs. 4.7(c) and 4.7(d) respectively. The results for scenario B are shown in Fig. 
for scenario C are shown in Fig. 14.91 

In order to quantify the behaviour, I give the values for the signal and background cross 
sections, the significance S and the signal to background ratio r for a specific set of beam polar- 



izations {Pe+\Pe-) = (0|0), (0|0.8), (-0.3|0.8), (-0.6|0.8), (0|0.9), and (-0.3|0.9) inTables[ 
I4.5[ 14.61 for the scenarios A, B, and C, respectively. I find that an additional positron polariza- 
tion Pg+ = —30% enhances the significance S by factors {1.5,1.5,1.6} in scenarios {A,B,C}, 



respectively, compared to beams with only e polarization {Pe+\Pe-) = 
tors {1.4,1.5,1.5} in scenarios {A,B,C}, respectively, for {Pe+\Pe-) = (■ 



0|0.8), and by fac- 
-0.3|0.9) compared 



to (P, 



(0|0.9). The signal to background ratio r is enhanced by {1.7,1.7,1.8} for 
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Table 4.4.: Cross sections a, significance S, and signal to background ratio r for different beam 
polarizations {Pe-\Pe+) for Scenario A at a/s = 220 GeV, with C = 500fb"\ 



Scenario A 


(0|0) 


(0|0.8) 


(-0.3|0.8) 


(-0.6|0.8) 


(0|0.9) 


(-0.3|0.9) 




6.7fb 


12 fb 


16 fb 


19 fb 


13 fb 


16 fb 


cr{e^e^ vvj) 


2685 fb 


652 fb 


534 fb 


416 fb 


398 fb 


360 fb 


S 


2.9 


10 


15 


20 


14 


19 


r 


0.3% 


1.8% 


2.9% 


4.6% 


3.2% 


4.6% 



Table 4.5.: Cross sections a, significance S, and signal to background ratio r for different beam 
polarizations {Pe-\Pe+) for Scenario B at a/s = 350 GeV, with C = 500 fb^^. 



Scenario B 


(0|0) 


(0|0.8) 


(-0.3|0.8) 


(-0.6|0.8) 


(0|0.9) 


(-0.3|0.9) 


(T{e+e- XiXil) 


5.5 fb 


9.6fb 


13 fb 


15 fb 


10.2 fb 


13.3 fb 


a{e^e^ vvj) 


3064 fb 


651 fb 


481 fb 


312 fb 


350 fb 


272 fb 


S 


2.2 


8.4 


13 


19 


12 


18 


r 


0.2% 


1.5% 


2.6% 


4.9% 


2.9% 


4.9% 



Table 4.6.: Cross sections a, significance S, and signal to background ratio r for different beam 
polarizations {Pe-\Pe+) for Scenario C at ^/s = 500 GeV, with C = 500 fb^^. 



Scenario C 


(0|0) 


(0|0.8) 


(-0.3|0.8) 


(-0.6|0.8) 


(0|0.9) 


(-0.3|0.9) 


cr(e+e- ^ XiXi^) 


4.7fb 


8.2 fb 


11 fb 


13 fb 


8.6 fb 


11.2 fb 


a{e^e^ vvj) 


3354 fb 


689 fb 


495 fb 


301 fb 


356 fb 


263 fb 


S 


1.8 


7 


11 


17 


10 


15 


r 


0.1% 


1.2% 


2.2% 


4.4% 


2.4% 


4.3% 



{Pg+\P,-) = (-0.3|0.8) compared to {Pe+\Pe-) = (0|0.8) and by {1.4,1.7,1.8} for (Pe+|P,-) = 
(— 0.3|0.9) compared to {Pe+\Pe-) = (0|0.9). If the positron beams would be polarized by 
Pg+ = —60%, the enhancement factors for S are {2,2.3,2.4}, and for r they are {2.5,3.2,3.6}. 
For Pg- = 0.8, it is only with positron polarization that I obtain values of r clearly above 1%. If I 
have Pg- = 0.9, then r exceeds 1% without positron beam polarization. 

Since the neutralinos are mainly bino, the signal cross section also depends sensitively on the 
mass mgj^ of the right selectron. In scenarios {A, B, C} the masses are = {133, 191, 261} GeV, 
respectively, see Table 14.31 For larger masses, the signal to background ratio drops below 
r < 1%. With {Pe+\Pe-) = (-0.3|0.8), this happens for m,-^ = {214,300,390} GeV, and the 
significance would be S < 5. These selectron masses correspond to the mSUGRA parameter 
Mo = {190,270,350} GeV. 
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4.4.3. Summary and Conclusions 

I have studied radiative neutralino production e+e^ — > XiXil the ILC with longitudinally 
polarized beams. For the center-of-mass energies ^/s = 220 GeV, 350 GeV, and 500 GeV, I have 
considered three specific mSUGRA inspired scenarios. In my scenarios, only radiative neu- 
tralino production is kinematically accessible, since the other supersymmetric particles are too 
heavy to be pair produced. I have investigated the beam polarization dependence of the cross 
section from radiative neutralino production and the background form radiative neutrino pro- 
duction e+e^ vv^. 

I have shown that polarized beams enhance the signal and suppress the background simul- 
taneously and significantly. In my scenarios, the signal cross section for (Pg+ \Pe-) = (— 0.3|0.8) 
is larger than 10 fb, the significance S > 10, and the signal to background ratio is about 2 — 3%. 
The background cross section can be reduced to 500 fb. Increasing the positron beam polariza- 
tion to Pg+ = —0.6, both the signal cross section and the significance increase by about 25%, in 
my scenarios. For (Pg+ |Pe-) = (0.0|0.9) the radiative neutralino production signature is observ- 
able at the ILC but both the significance and the signal to background ratio are considerable 
improved for (Pe+lPg ) = (— 0.3|0.9), making more detailed investigations possible. The elec- 
tron and positron beam polarization at the ILC are thus essential tools to observe radiative 
neutralino production. For unpolarized beams this process cannot be measured. 

I conclude that radiative neutralino production can and should be studied at i/s = 500 GeV, 
as well as at the lower energies ^/s = 220 GeV and ^/s = 350 GeV, which are relevant for Higgs 
and top physics. I have shown that for these energies there are scenarios, where other SUSY 
particles like heavier neutralinos, charginos and sleptons are too heavy to be pair produced. 
In any case, a pair of radiatively produced neutralinos is the lightest accessible state of SUSY 
particles to be produced at the linear collider. 



4.5. Summary and Conclusions 

I have studied radiative neutralino production e+e^ A^i^iT at the linear collider with po- 
larised beams. I have considered the Standard Model background process e+e^ 1/V7 and the 
SUSY background e+e^ — > 1/17*7, which also has the signature of a high energetic photon and 
missing energy, if the sneutrinos decay invisibly. For these processes I have given the complete 
tree-level amplitudes and the full squared matrix elements including longitudinal polarisations 
from the electron and positron beam. In the MSSM, I have studied the dependence of the cross 
sections on the beam polarisations, on the gaugino and higgsino mass parameters M2 and }i, as 
well as the dependence on the selectron masses. Finally, in order to quantify whether an excess 
of signal photons, Ns, can be measured over the background photons, Nb, from radiative neu- 
trino production, I have analysed the theoretical statistical significance S = Ng/ yjNs + Nb and 
the signal to backgroimd ratio Ns/N^- Our results can be summarised as follows. 

• The cross section for e+e^ XiXil reaches up to 100 fb in the ii-M2 and the mo-M2 plane 
at ^/s = 500 GeV. The significance can be as large as 120, for a luminosity of £ = 500 fb^ , 
such that radiative neutralino production should be well accessible at the ILC. 
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• At the ILC, electron and positron beam polarisations can be used to significantly enhance 
the signal and suppress the background simultaneously. I have shown that the signif- 
icance can then be increased almost by an order of magnitude, e.g., with {Pe-,Pe+) — 
(0.8,-0.6) compared to (Pe-,Pg+) = (0,0). In the SPS la scenario the cross section C7"(e+e^ - 
/^i/ti7) increases from 25 fb to 70 fb with polarised beams, whereas the background (i{e^ 
1/V7) is reduced from 3600 fb to 330 fb. Although a polarised positron beam is not essential 
to study radiative neutralino production at the ILC, it will help to increase statistics. 

• I note that charginos and heavier neutralinos could be too heavy to be pair-produced 
at the ILC in the first stage at a/s = 500 GeV. If only slepton pairs are accessible, the 
radiative production of the lightest neutralino might be the only SUSY process to study 
the neutralino sector. Even in the regions of the parameter space near the kinematical 
limits of - P^i^ production I find a cross section of about 20 fb and corresponding 
significances up to 20. 

• Finally I want to remark that my given values for the statistical significance S can only 
be seen as rough estimates, since I do not include a detector simulation. However, since 
I have obtained large values up to S ~ 120, I hope that my results encourage further 
experimental studies, including detailed Monte Carlo simulations. 
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Figure 4.7.: Signal cross section (a), background cross section (b), significance (c), and signal to 
background ratio (d) for a/s = 220 GeV, and an integrated luminosity C = 500 fb^^ 
for scenario A: Mq = 90GeV, M1/2 = 225 GeV, Aq = -90GeV, and tanjS = 10, see 
Tables and m 
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Figure 4.8.: Signal cross section (a), background cross section (b), significance (c), and signal to 
background ratio (d) for a/s = 350 GeV, and an integrated luminosity C = 500 fb^^ 
for scenario B: Mq = 135 GeV, M1/2 = 325 GeV, Aq = -135 GeV, and tanjS = 10, 
see Tables i^l and m 
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Figure 4.9.: Signal cross section (a), background cross section (b), significance (c), and signal to 
background ratio (d) for a/s = 500 GeV, and an integrated luminosity/! = 500 fb^^ 
for scenario C: Mq = 200 GeV, M1/2 = 415 GeV, Aq = -200 GeV, and tanjS = 10, 
see Tables 1121 and Us] 
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5.1. Introduction 

If supersymmetric particles are discovered, the underlying SUSY parameters can be determined 
from measurements of cross sections, particle masses, decay widths, and branching ratios. 
Many authors developed methods and programs to extract the parameters from these mea- 
surements. In the following, I shall give an overview over the methods, concentrating on the 
gaugino sector. 

Choi et al. analyse in [109] the chargino system. They present an analytical method to extract 
the parameters M2, ji, and tan /3 of the chargino mixing matrix from chargino pair production 
in e+e^ annihilation with polarized beams. The absolute errors on M2 and j^ are of the order of 
GeV, and the error on tan /3 is if tan j6 is not too large. 

In [40] the analysis has been extended to the neutralino system to obtain the bino mass param- 
eter Ml . As I shall demonstrate later, this method demands chargino parameters and neutralino 
masses measured with an accuracy 0(0.1 GeV), which is not feasible in the first run of the ILC. 

Desch et al. present in [110-112] a study to determine the parameters Mi, M2, w, and tanj6 
from a fit of the light and heavy neutralino and chargino masses to LHC and LCU data. They 
present formulae to determine M2, ff, and tan j6 from the chargino masses and from cross sec- 
tions with left (P+|P-) = (0.6| -0.8) and right (P+|P-) = (-0.6|0.8) longitudinally polarised 
beams. Mi is obtained from the polarised cross sections cr(e+e^ — > X1X2) ^{e^e^ — > 
;^2^2)- They simulated the parameter determination with an LC measurement at the SPSla 
point [83,84]. They could recover the input data with absolute errors 0(0.1 GeV) for Mi and 
M2, 0(1 — lOGeV) for }i, and 0(1) for tanj6. Combining the analysis with LHC data reduces 
the errors on these parameters by a factor of about 2. 

Bechtle et al. present in [113] the program Fittino. It performs non-linear fits to observables 
such as masses, cross sections, branching fractions, widths and edges in mass spectra to deter- 
mine the SUSY parameters. More about non-linear fits can be found in Ref. [114]. The authors 
implemented an iterative fitting technique and the simulated annealing algorithm to obtain the 
fit parameters. In [112] they present an example calculation using the SPSla point. They re- 
covered the input data with errors of about 0(0.01 GeV) for Mi and M2, 0(1 GeV) for f/, and 
O(O.l) for tan|6. 

In [115] the authors discuss the parameter determination in a focus point inspired scenario. 
The slepton and squark masses are about 2000 GeV, which is even heavier than the particle 
spectrum of the SPS2 point. The determination of Mi and M2 succeeds with an error of 0(0.1 — 
1 GeV), but the errors on }i and tan /3 are 0(10 — 100 GeV) and 0(10), respectively. 

Sf itter [116] is another program to extract SUSY parameters from particle masses using a 
fit or multi-dimensional grid or both. The fit is performed assuming the mSUGRA parameters. 

^In 2002, the terminology was "LC". Today, we are talking about the "ILC". 
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In a simulation with the SPSla point they get errors of about 0{1 GeV) for Mi, M2, and ^; the 
error on tan j6 is ~ 3. 

In [117] the authors use Markov chain techniques to determine SUSY mass measurements 
from simulated ATLAS data. 

The Supersymmetry Parameter Anaylsis project (SPA) [36] provides a common framework 
for parameter determination. The authors present a list of computational tools to perform the 
required calculation: These are tools to translate between calculational schemes, spectrum cal- 
culators, calculators for cross sections, decay widths etc., and event generators, parameter anal- 
ysis programs, RGE programs, and auxiliary programs. The authors define the tasks of the SPA 
project as follows: promoting higher order SUSY calculation, improving the understanding of 
the DR scheme, improving experimental and theoretical precision, improving coherent analyses 
from LHC and future ILC data and determining SUSY parameters, determining and clarifying 
the nature of dark matter, and the study of extended SUSY scenarios. 

AUanach et al. [118] use genetic algorithms to distinguish between different SUSY models. 

All the methods above determine mass parameters and fundamental SUSY parameters such 
as Ml, M2, and tan j5. The errors on the gaugino parameters are small, the errors on the 
higgsino mass parameter and on tan j6 are somewhat larger. 

I present a method that determines the couplings of the lightest neutralino. From theses cou- 
plings, I calculate the corresponding elements of the neutralino diagonalisation matrix, assum- 
ing unitarity. The absolute errors on the elements of the diagonalisation matrix are of order 
0.001 — 0.01. With the knowledge of the neutralino masses, I then obtain the values of Mi, 
M2, ji, and tanj6. The errors are 0.4 GeV, 4 GeV, 2.5 GeV, and 7, respectively. This method is 
complementary to the methods described above. This allows for cross checks. 



5.2. The circle method 

The authors of [40] present a method to calculate Mi and (pM^ for the CP violating extension of 
the MSSM from the characteristic polynomial of the matrix il.2} . Mi and }i are here complex 
parameters, cf Eq. (|1.4|) : 

= det(M+M - m|o) = fn|o - flwjjo + fcwjjo - cmj.o + d, (5.1) 

with the polynomial coefficients given by 

a = \Mi\^ + Ml + l\}i\^ + lml, (5.2) 
b = |MipM| + 2|/^|2(|Mi|2 + M^) + (|/^P + m|)2 
+2?n| { |Mi 1^ cos^ e^o + Ml sir? e^u- 

sin2,6 [|Mi I sin^ Ou, cos(<^i + ^2 cos^ 0,^ cost/),,] } , (5.3) 

c = \^f {\^\'^{\Ml\^ + Ml)+2\M^fMl + mzsm^2^ + 2mz{\Ml\^cos^eu, + Mlsm^9-a,)} 

- 2m| I I sin 2,6 { I Ml I (M^ + I 1 2 ) sin^ cos ((^1 + ) + M2 ( I Ml I ^ + I I ^ cos^ ) cos } 

+m| I |Mi p cos^ 6,, + 2|Mi IM2 sin^ 6,, cos^ 6,, cos cpi + Mf sin^ 9,^ } , (5.4) 
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d = I Ml pMf I I ^ - 2m|| p I Ml I M2 sin 2^ { I Ml 1 2 cos^ cos + M2 sin^ 9,, cos (^1 + } 

+m||/ipsin2/3||Mipcos'^0i„ +2|Mi|M2sin2 0i(,cos^0i(,cos(^i + M2sin"^0a,| . (5.5) 

Note that the matrix M is symmetric but not hermitian. So one has to diagonaUse M+M to 
get the singular values of M which are the physical masses. Eq. i5.1} is quadratic in 9^eMi 
and DimMi for fixed m^o, i = \ . . A, and for fixed parameters M2, }i, and tan /3. So it describes 
four circles in the D^eMi - JmMi -plane, which should intersect in one point (D^eMi, 3mMi). In 
general, four circles do not intersect in only one point. There may be no solution or even two 
solutions, in which case all midpoints have to be located on a straight line. From the coordinates 
of this point one can calculate |Mi| and (pu^ = arg(Mi): 

= det(M+M-m^o) 

= A(m^o,M2,f/,tan|6)X^ + A{m^o, M2, pi,tan ^)Y^ 

+Bi (m^o, M2, tan j6)X + B2{m^o, M2, f/, tan j6)y - C (ra^o, M2, pi, tan j6) (5.6) 

with X = ^HeMi, Y = 3mMi. The radii r/ and the midpoints m, of the circles described by 
Eq. (|5.6|) are given by 



mi{m,\my) = M,- (^-— -—J 



(5.8) 



In Fig. 5.1(a) I show the four circles for all neutralino mass m^o. As input data I have taken 



the RP" model of Ref. [40]: 

(|Mi|,(/>i,M2,f/,(/?,„tan|6) = (100.5 GeV, f, 190.8 GeV, 365.1 GeV, f, lO) , (5.9) 
leading to the following neutralino masses 

(m^o,m^o,m^o,m^o^ = (99.15 GeV, 177.07GeV, 372.0 GeV, 387.41 GeV) . (5.10) 

It is clear that the four circles intersect in one point: X = 52.8 GeV, Y = 85.5 GeV. This yields 
|Mi| = 100.5 GeV and cpi = 1.02 7r/3. These two values agree with the input data (|5.9|) . 
The algebraic form of the coefficients A, Bi, B2, and C follow from the Eqs ( I5.1| )-( I5.5D . For their 
numerical values one needs the values of the parameters M2, tan /5, and (pfi, which can be 
determined from the chargino system, see Ref [40, 109], and at least three neutralino masses 
form LHC/ILC measurements. 

Measurements of masses and cross sections have unavoidable errors. These errors influence 
the radii and the midpoints of the circles Eq. (|5.6|) . I have analysed how errors in the parame- 
ters M2, tan/3, and (pfi and the neutralino masses m^o influence the circles of the example 

given in [40] and found that the circles belonging to the neutralinos X2 ' xt drift away con- 
siderably from their exact position even for small errors. This behaviour is due to zeros and 
poles of the coefficients a = a{M2,}i,tanfi), hi = bi(M2, /i,tan/3), ^2 = f'2(M2, f/, tan/3), and 
c = c(M2, tan/3) which determine the radius and the midpoint of the neutralino circles. In 
Fig. 5.2(a)|(572(d) I illustrate this behaviour of the radius for the neutralino mass circles of xt- 
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-200 200 400 

(a) Exact circles 




-200 4-r 



-200 



200 400 



(b) Circles with = -t-0.5GeV 



Figure 5.1.: Influence of a small error in ^ on the intersection point of the four circles in the 
9^eMi - 3mMi - plane. The colors mean: red: cyan: x.2' blue: x^' green: x^- 



The poles and the zeros are located very close to each other and close to the neutralino masses. 
This leads to the messy situation that small errors in the input data blow up to large errors in the 
radii and coordinates of the midpoints. The neutralino mass matrix (|1.2|) can be decomposed in 
main diagonal and off -diagonal blocks. The off -diagonal blocks are proportional to mz- Then, at 
zeroth order in niz, the eigenvalues of M are given by m^o ~ |Mi 



m 



X2 



M2, and m^o 

A3/4 



these eigenvalues are not necessarily mass ordered. This relation helps to understand why the 
mass circle of the fourth neutralino reacts so strongly on the errors on m^o and }i. 



In Fig. 5.1(b) I show for the RP" model [40] how the circles drift away if the measured value 
of ji is 0.5 GeV larger than the "true" value. This corresponds to a 0.1% error! There is neither 
an intersection point nor a small region where all possible pairs of neutralino circles intersect. 

In Fig. l5.3l and Fig. l5.4[ 1 show how the circle of the fourth neutralino is disturbed by errors in 
the values of m^o, M2, }i, tan/5, and (pji. One of these five parameters is varied, the others are 
kept at their exact values. 

In each picture of Fig. 15.31 and Fig. I5.4[ I show three circles: two perturbed circles and one 
unperturbed (black) circle. In one case the parameter is a little bit too large (green circle) , in the 
other case a little bit too low (red circle). The other parameters are not varied. The unperturbed 
circle (black) is thus the same in all Figs (note however the scale change). 

From these figures 1 conclude for the circle of the fourth neutralino: 

• ra^o must be measured very precisely, small errors lead to large deviations from the un- 
perturbed circle. The upper bound on this error is Ara^o ~ 0.1 GeV. 

• The dependence on M2 is weaker. A maximal error on about AM2 ~ 1 GeV is allowed. 
The reason for the weak dependence is that m^o does not depend on M2 to zeroth order. 

• The experimental error on }i should be small: Aji < 0.1 GeV. This is due to the strong 
dependence of m^o on ji. For the error on its phase, I find: Acp^ < 0.01 tt. 
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800- 




OH — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 
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(a) Radius as a function of ^ for the correspond- 
ing neutralino mass circle. The exact value for ji is 
365.1 GeV, and the mass of xl is w^o = 387.4 GeV, see 
Eqs and <5.10t . The other parameters are kept on 
their exact values. 




170 180 190 200 210 220 

Ml peV] 



(b) Radius r4 as a function of M2 for the correspond- 
ing neutralino mass circle. The exact value for M2 is 
190.8 GeV, and the mass of xl is m^o = 387.4 GeV, see 
Eqs J5.91 and JS.lOt . The other parameters are kept on 
their exact values. 
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800- 




360 370 380 390 400 410 
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(c) Radius as a function of m^o for the correspond- 
ing neutralino mass circle. The mass of X4 is fn^o = 
387.4 GeV, see Eqs (SJ) and iSlOt . The other parame- 
ters are kept on their exact values. 




40 80 120 160 



(d) Radius as a function of (pp for the corresponding 
neutralino mass circle. The exact value for ip is 7t/4, see 
Eqs l l5.9t . The other parameters are kept on their exact 
values. 



Figure 5.2.: Dependence of the radius on ji, M2, m^o, and (p^ 
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/1_ \ AST n 1 _"\ 7 



n 1 _^ ?■ 









(k) Am^o = 1.0 GeV 



(h) AM2 = 1.0 GeV 



(i) A?( = 1.0 GeV 






(j) Am^o = 2.0 GeV 



(k) AM2 = 2.0GeV 



(1) Afi = 2.0 GeV 



ure 5.3.: Influence of small perturbations of m^o (left column), M2 (middle column), and on ji (right 
column) on the circle of neutralino Xir black: unperturbed circle, green: +A[m^o,M2,fi), 

A4 



red: — A()«-o,M2, fO- black-white-printing, black = black, red = dark-grey, green = light 



grey. The color code is for all figures the same as in Fig. 5.3(i) 
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(g) A(pf, = 0.1 TT (h) Atan^ « 3.17 

Figure 5.4.: These figures show how errors on (ppi and tan /3 influence the position of the circle of 
neutralino ;f 4. black: unperturbed circle, green: +A(<^^„tanj6), red: — A(^^,, tanj6). 
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• For the error on tan j6, 1 find A tan /3 < 0.3. 

• The expected experimental errors of M2 and }i at the ILC are 0{1 GeV), see Ref [109]. This 
error is too large for the method described in Ref [40]. 

The situation is similar or worse for the circles of the neutralinos Xi Xy For the circles 
of the second neutralino the influence of the error on M2 is disastrous. I conclude that it is not 
possible to determine |Mi | and (pi with the circle method from Ref [40]. In the following section 
I propose an alternative method. 

5.3. Determining Neutralino Couplings 

I will show in this sections how radiative neutralino production together with neutralino pair 
production 

e^e- ^ xU' (5-11) 

can be used to determine the couplings of the neutralinos to the Z boson and to Cr/l- My 
method is not only designed for the MSSM, it is applicable to every model with measurable 
cross sections. The idea is as follows: Write the cross section a of an arbitrary process as 

a = Yici{ai,bi,fi)Xi , (5.12) 

i 

where the c, are functions of the unknown couplings fl;, and /y, the functions Cj are not 
necessarily linear. The X, are calculable factors, depending only on the neutralino and selectron 
masses. If there are n cross section measurements, n > number of couplings, then one can 
perform a least square fit to determine the couplings fl;, bi, and /y . If the Eq. ([5.12|) is non- 
linear in the couplings one can either linearize Eq. ([5.12|) and use the custom linear least square 
functions provided by Maple or Mathematica or one can use techniques for non-linear fits like 
Minuit [119]. I have chosen the first approach because the linearized equations are not too 
complicated. The method is best illustrated by an example. 

5.3.1. Mathematical Structure of the cross section and the couplings 

To determine the couplings of the neutralinos to the selectrons and the Z^-boson, I use the cross 
sections of the following reactions: 

e+e- ^ xUh (5-13a) 
e^e- ^Xixl z = 2,3,4 , (5.13b) 
e^e- ^ xlxl ■ (5.13c) 

It is straightforward to include further reactions e+e^ — > X^X^, ij = 23,24,33,34,44, if they are 
measurable. Their cross sections can be decomposed as follows: 

a{e+e- ^ xlxh) = a^^ = a\X + h\Y + FiZ, (5.14a) 

a{e+e- ^ xW) = = alajXy + blbjYu + flifl,/i,X2,- - foifo./i^Yz; + /i';Z;, (5.14b) 

(T{e+e- ^ X2X2) = ^21 = 4^1,22 + fo!ri,22 + fli/22X2,22 " ^1/22^2,22 + /2'2Z22 • (5.14c) 
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The factors fl,, are associated with the couplings to the right-selectron, left-selectron, and 
Z°-boson, respectively, and are given by 

(5.15a) 



cos 010 ' 



1 / N,2 N,i , 



2 V sin 0„, cos I 

/y = ^ (N^Np - N,-4N^-4) , (5.15c) 
Fi=/i'i. (5.15d) 

/ii appears only quadratically in Eq. (|5.14a|) and may be small, its error from a least square fit 
however may be large. The optimal value for /^^ can become negative, which is unphysical. 
Therefore, the use of Fi instead of /n as a fit parameter secures the convergence of the iteration. 
Due to this problem, the value of f\\ is not used further. The X„ Y,, Z, are functions of the 
right selectron mass, left selectron mass, and Z'^-mass, respectively. They all depend on the 
neutralino masses, the center of mass energy and the longitudinal polarisation of the electron- 
positron-beam. Their explicit form can be found in [120]. They need to be calculated only 
once. So one does not have the problem that the iterations in the program that tries to find the 
minimum of 'X^ does not converge due to errors occurring when integrating out the X,, Y,, and 
Z; by Monte Carlo integration. 

This set of equations is nonlinear in the couplings parameters. The equations can be expanded 
in a Taylor series up to first order: 

'^(P) ~ ^(Po) +t^(po)(p- Po). (5.16) 

where p is a vector, collecting the parameters fl,, hi, fjj. The linearized equations can be solved 
by a least square fit recursively, po is a first guess of the solution. 

From the parameters a,, bi, and fij one can determine the entries of the neutralino diagonali- 
sation matrix N. 

Nil = - cos 0i„fl,, 2 = 1 ... 4, (5.17a) 
Ni2 = sin {2bi + ai), i = \..A, (5.17b) 



N23 = ±\l\{l- N|i - + 2/22), (5.17c) 



N24 = ±^^(1 - - - 2/22),= - Nl, - - Nfg), (5.17d) 

2/12 - Nil N21 - N12N22 ,^-,7, 

= ms ' ^'-''"^ 

2/12 + Nil N21 + N12N22 

Ni4 = 2N2^ , (5.17f) 

2/13 - N11N3I - N12N32 

^33 = 2Ni^ . (5.17g) 
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N34 = ± \/ 1 - Nf, - - NI3 , (5.17h) 



2/14 - N11N41 - N12N42 ,^ , 

W43 = ^TT , (5.1/1) 



N44 = ± ^ 1 - N|i - - NI3 , (5.17j) 

Eqs. I|5.17a|l - (|5.17b|l are derived from Eqs. (|5.15a|) - (|5.15b|l , N23 and A/24 are obtained from the 
unitarity relation of the matrix N; f22, N13 and N14 are constructed in such a way, that = 
(Nii,Ni2, Ni3, N14) is orthogonal to ^ = (N21, N22, N23, N24); the elements N33, N43 are cal- 
culated from /i3 and fu as well as the orthogonality relation ■ ^ = 0. We use unitarity 
to calculate N43 and A/44, because this leads to a smaller error for these elements. The sign is 
choosen in such a way, that ^ are orthogonal to Xi- 



5.3.2. The cross sections for Neutralino pair production 

In Fig. I5.5[ I show the cross sections for e+e^ — > XiX^^ ^ = 2 ... 4 and e+e^ XiXil for 
cms - energies from 200GeV - lOOOGeV for three different polarisations: (P+|P_) = (0|0) in 
Fig. [53(a)| (P+|P-) = (-0.6|0.8) in Fig. [53(b)l and (P+|P-) = (0.6| -0.8) in Fig. [53(c)| The 



figures show, how suitable beam polarisation enhances cross sections. (P_|_|P_) = (0.6| — 0.8) 
enhances x^Xi' X2X2' P^i^ production compared to unpolarised beams, the opposite beam 
polarisation enhances radiative neutralino production, X1X2 and X1X3 and XiXt- Iri my example 
model which I will present below in detail, the Xi is mainly bino (~ 95%) which couples mostly 
to right handed sleptons, the X2 mainly wino (~ 85%) which couples preferably to left handed 
slepton; so cross sections with Xi involved are enhanced by right handed beam polarisation, 
and cross sections with X2 are enhanced by left handed beam polarisation. 

Polarised beams are essential for the described method to determine parameters, since they 
enhance couplings either between right handed particles or left handed particles. 



5.4. An example 
5.4.1. The model 

In order to demonstrate my method, I choose an example model with light neutralinos, so that at 
least the processes e+e^ xlXil' ^^^^ ~^ XiX^' ^ = 2,3,4, and e+e^ X2X2 kinematically 
accessible at the ILC and the cross sections exceed 0(10 fb) for both polarisations. The higgsino 
components of Xi should not be too small, so that there might be a chance to determine the 
Xi-Xi-Z- coupling. 

The input data and the derived neutralino and selectron masses are listed in the first row 
of Tab. 15.11 For comparison I also list the values for the SPSla scenario. In my model the cross 
section for the process e^e^ — > X2X3 is larger than 10 fb for both beam polarisations and the cross 
section of the process e^e^ — > X2X4 is larger than 10 fb for left beam polarisation. Therefore, I 
shall later present a study including these processes. The fairly light particle spectrum and, 
additionally, the large wino component in the latter process are the reasons for the large cross 
sections. I do not include them from the beginning. 
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(c) Beam polarisation (P+|P-) = (0.6| - 0.8) 
Figure 5.5.: Comparison of cross sections for different beam polarisations. 
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M2 


Ml 




Mo 


tan ,6 




A2 


m-.o 

A3 








my model 


165 


82.5 


230 


75 


11 


76.3 


136.8 


-239.2 


273.2 


117.3 


171.3 


SPSla 


192 


102 


352 


100 


10 


99 


175 


348 


369 


145 


204 



Table 5.1.: Input and mass parameters of the example model. The minus-sign appearing in the 
third neutralino mass denotes the CP-eigenvalue of this particle. All masses are given 
in GeV. 



With these values the neutralino diagonalisation matrix follows as 



N = 



f 


0.953 


-0.117 


0.257 


-0.106\ 




-0.241 


-0.844 


0.402 


-0.262 




-0.089 


0.129 


0.682 


0.715 




-0.158 


0.508 


0.554 


-0.640/ 



(5.18) 



The neutralino couplings to eR,ei, and Z° are listed in Tab. 



couplings to Sr couplings to ei couplings to 



«1 = 


-0.740 




0.287 


/ii = 


0.065 


fl2 = 


0.187 


b2 = 


-0.690 


/12 = 


0.038 


«3 = 


0.069 


b3 = 


0.057 


/l3 = 


0.126 


fl4 = 


0.123 


h = 


0.398 


fu = 


0.037 










fll = 


0.047 



Table 5.2.: Theoretical couplings. 



I assume that eight cross section measurements of each process are available: Four different 
cms energies of the beam (500 GeV, 550 GeV, 600 GeV, 650 GeV) are combined with two different 
longitudinal beam polarisations (P+|P-) = (-0.6| +0.8) and (P+|P-) = (+0.6| -0.8). 

Each measurement has an error. The error on the cross sections consists of the statistical Pois- 
son error and the systematic error. I consider only the statistical error. To simulate the statistical 
error on a measurement, 1 calculate the exact cross section and add a Gaussian distributed ran- 
dom number with zero mean and variance V = (Sa)'^. The statistical error Scr follows from 

^e+e-^xUh = ^e+e-^By " K+e-^vvj' ^process = Cr{pmceSs)C, 
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~'Ai A j 



(5.20) 



with C and Nprocess denoting the integrated luminosity in fb^^ and the number of events of 
the process, respectively. Eq. \5.19) accounts for the fact that the number of events for radiative 
neutralino production is calculated as the difference from all events "photon plus missing en- 
ergy" and the radiative neutrino background. Eq. (|5.20|) is the Poisson error for neutralino pair 
production. 

The values for the cross sections and their errors are: 





Polarisation 




ai2 


C^13 


(7i4 


0-22 


[GeV] 


in 


% 




all cross sections in fb 




500 


(-0.6| 


+ 0.8) 


75.6 ± 1.0 


118.3 + 0.5 


81.0 + 0.4 


43.4 + 0.3 


10.5 + 0.1 


500 


(+0.6| 


-0.8) 


(3.4 + 4.5) 


181.5 + 0.6 


10.0 + 0.1 


17.2 + 0.2 


328.7 + 0.8 


550 


(-0.6| 


+ 0.8) 


68.3 ± 1.0 


105.6 + 0.5 


67.2 + 0.4 


42.6 + 0.3 


10.1 + 0.1 


550 


(+0.6| 


-0.8) 


(3.11 + 4.6) 


166.5 + 0.6 


8.6 + 0.1 


17.7 + 0.2 


317.9 + 0.8 


600 


(-0.6| 


+ 0.8) 


61.7 + 1.0 


94.0 + 0.4 


56.2 + 0.3 


40.0 + 0.3 


9.4 + 0.1 


600 


(+0.6| 


-0.8) 


(2.8 + 4.6) 


151.7 + 0.5 


7.5 + 0.1 


17.2 + 0.2 


300.4 + 0.8 


650 


(-0.6| 


+ 0.8) 


55.8 + 0.9 


83.9 + 0.4 


47.5 + 0.3 


36.8 + 0.3 


8.8 + 0.1 


650 


(+0.6| 


-0.8) 


(2.6 + 4.6) 


137.8 + 0.5 


6.6 + 0.1 


16.2 + 0.2 


280.3 + 0.8 



(5.21) 



The values in brackets are not used for further calculations as they are too small. These input 
data from Tab. 15.211 lead to 36 equations for 13 fit parameters. The diagonalisation matrix has 
16 entries, but only six of them are independent because of unitarity. There are ten relations 
between the matrix elements together with 13 equations from the fit parameters. The system is 
overdetermined. I choose the equations to determine the elements of the diagonalisation matrix 
such, that the error on the elements is as small as possible. In principle, unitarity can be tested 
by the first and the second column. So it can be tested if there is an additional singlino field. 

For the couplings I get as a result of the least square fit: 



couplings to 



couplings to 



couplings to 



fll = 


-0.7397 + 0.007 


h = 


0.2882 + 0.003 


{fk = 


-0.035) 


fl2 = 


0.1879 + 0.003 


b2 = 


-0.6918 + 0.002 


hi- 


0.036 + 0.006 


«3 = 


0.0667 + 0.001 


h = 


0.064 + 0.002 


/13 = 


0.131 + 0.002 


fl4 = 


0.122 + 0.002 




0.303 + 0.006 


/14 = 


0.043 + 0.002 










fll = 


0.051 + 0.009 



Table 5.3.: Result of the fit on the couplings 



The obtained ^ value is 



20.9, 



X^ 
36-13 



0.91 



(5.22) 
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With these data I get as a neutralino mixing matrix: 



/ 



N 



0.953 
-0.242 
-0.086 
-0.157 



-0.116 
-0.845 
0.138 
0.514 



0.252 
0.405 
0.714 
0.585 



-0.124 \ 
-0.249 
0.681 
-0.607 / 



± 



/ 0.009 0.005 0.01 0.05 \ 

0.004 0.002 0.015 0.04 

0.001 0.003 0.030 0.031 

\ 0.002 0.008 0.035 0.030 / 



(5.23) 



The error matrix shows that the bino and the wino components can be determined with high 
accuracy. The higgsino components have large(r) errors, especially the ones of the heavy neu- 
tralinos. 

The parameters Mi, M2, and }i follow from 



M = diag{mi,m2,m3,m^o)N 



Ml = Mil = 82.8 ± 1.2 GeV, 
M2 = M22 = 166.5 ± 2.3 GeV, 
ji = -M34 = 229.6 ± 2.5 GeV 



82.8 


0.38 


-5.35 


39.2 


0.38 


166.5 


9.4 


-77.8 


-5.35 


9.43 


-1.06 


-229.6 


39.2 


-77.8 


-229.6 


-0.43 



(5.24a) 

(5.24b) 
(5.24c) 
(5.24d) 



I do not derive any value for tan /3 because the result is not very reliable. The relative error on 
Ni3, Ni4, N23, and N24 are about 10% (error propagation in this 2x2 sub-block), and together 
with the bad behaviour of the tan function near its poles this leads to a result with a large error. 

5.4.2. How much does radiative neutralino production improve the measurements? 

If I omit the data from the process e+e^ X\Xilf ^rid repeat the calculation, then I get for the 
couplings the following result: 



couplings to Sr 



couplings to Si 



couplings to 



fll = 


-0.736 ± 0.039 


h = 


0.278 ± 0.005 






fl2 = 


0.180 ± 0.009 




-0.694 ± 0.005 


/l2 = 


0.056 ± 0.007 


«3 = 


0.058 ± 0.003 


h = 


0.099 ± 0.003 


/13 = 


0.154 ± 0.002 


fl4 = 


0.128 ± 0.007 


h = 


0.296 ± 0.008 


/14 = 


0.026 ± 0.005 










/22 = 


0.025 ± 0.017 



Table 5.4.: Result of the fit on the couplings without the data from radiative neutralino produc- 
tion. 



The errors on the elements of the neutralino mixing matrix are increased by about a factor of 
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2 — 5. This leads to the following values of the gaugino parameters 

Ml = 82.5 ± 6.3 GeV, (5.25a) 
M2 = 160.5 ± 3.1 GeV, (5.25b) 
}i = 220.6 ±6.2 GeV . (5.25c) 

The error on Mi is five times larger than the error of the case that includes radiative neutralino 
production. The error on M2 is enlarged only by a small amount and the error of }i is more than 
doubled. 

Including the cross section of radiative neutralino production leads to smaller errors on Mi, 
M2, and fi. So it is worthwhile to examine radiative neutralino production at a future linear 
collider. 



5.4.3. The effect of including the production of further neutralino pairs 

As I mentioned at the beginning of this section, the cross sections for 

e^e- ^ xlx3 (5-26) 
and e+e- ^ Xixl (5-27) 

can exceed 10 fb, see Fig. 15.61 These two processes introduce two further couplings, /as and /24, 
for their definition see Eq. (|5.15c|) . They are not necessary for calculating the elements of the 
neutralino diagonalisation matrix and I do not use them for further calculations. Nevertheless, 
they provide useful information about signs, see the next subsection. N33, N34, N43 and N44 can 
be expressed as 

2/23 - N21N31 - N22N32 ,^ _„ . 

N33 = TT / (5.28a) 

A/23 

2/23 + N21N31 + N22N32 ,^ , 

^ N^, ' ^ ^ ^ 

2/24 - N21N41 - N22N42 _„ , 

N43 = , (5.28c) 

A/23 

2/24 + N21N41 + N22N42 ,^ ^„ 

N44 = — • (5.28d) 

A/24 

The additional processes reduce the errors on a, and bj, i = 2, 3, 4. As a consequence the errors 
on Ml, Ml, and jx are reduced by 20%, respectively. 



5.4.4. Resolving Ambiguities 

The system of Eq. ( 15.131 ) has eight fix points. The starting point determines that fix point to 
which the iteration will converge. Some of these fix points do not fulfill the unitarity condition, 
so these points are to be discarded. The signs of Nn is choosen as +1. With this choice the signs 
of fl;, / > 1, are fixed, if the solution is physical. The other signs of Nij are choosen such that the 
eigenvectors are orthogonal to each other. 
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Ecms 

(a) Unpolarised cross section 





Ecms 



(b) Beam polarisation (P+ 1 P_) = (-0.6|0.8) 




300 400 



500 600 700 

Ecms 



800 900 1000 



(c) Beam polarisation (P+|P-) = (0.6| -0.8) 
Figure 5.6.: Comparison of cross sections for different beam polarisations. 
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5.4.5. Unitarity 

Throughout this method, I assumed unitarity, which means that the neutralino system is com- 
plete. This assumption can be dropped to test if there is an additional singlino field [121]. The 
eigenvector of the first and the second neutralino are suitable candidates to test unitarity. But a 
detailed analysis is beyond the scope of this study. 

5.4.6. Further Studies 

The presented method can be extended to 

• the MSSM with a CP violating gaugino sector, 

• NMSSM to test unitarity, 

• the chargino sector to determine the matrices. 

• The circle method could be used for precision measurements of the chargino parameters. 

5.5. Conclusion and Summary 

In this chapter I presented a method to determine the Xi'X.r^R/L arid Xi'Xr^^ couplings. 

• The method from [40] to determine Mi does not work because the circles are too sensitive 
to errors of the input data. 

• It is possible to determine the discussed couplings from the polarised cross sections of 
radiative neutralino production and neutralino pair production with errors 0(0.001 — 
0.01). The masses of the neutralinos and the selectrons must be known from LHC/ILC 
measurements . 

• From the couplings one can determine the neutralino diagonalisation matrix. The errors 
on the elements are about 0(0.001 - 0.01). 

• From the neutralino diagonalisation matrix and the neutralino masses one can determine 
the neutralino mass matrix. The errors of Mi, M2, and ji are about 0{1 GeV). It is difficult 
to determine tan fi with my method. 

• Omitting the cross sections of radiative neutralino productions enlarges the errors on Mi, 
M2, and Additional processes such as e+e^ X2X3/4: reduce the error on the corre- 
sponding couplings and on Mi, M2, and }i. 

• The differences to a global fit approach are as follows; In Fittino [113], the parameters 
Ml, M2, }i, and tan/3 are fitted directly to the data, and the couplings are obtained as a 
by-product. In Fittino, up to 24 SUSY parameters can be fitted to cross sections, edge po- 
sitions, branching fractions, cross sections times branching fractions and Standard Model 
parameters from LHC/ILC measurements. The authors recover Mi and M2 with absolute 
errors of the order 0(0.01 — 0.1 GeV), and }i with an absolute error of the order 0(1 GeV) 
(129 degrees of freedom). 

In my method, the couplings are the fit parameters to cross sections from ILC measure- 
ments. From the couplings, I obtain Mi, M2, and ji. The couplings enter the tree level 
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cross sections of the considered particles, which are in my case neutralinos. I recover the 
input parameters of my model with absolute errors of the order 0(1 GeV) (23 degrees of 
freedom). The coupling independent terms of the cross sections need to be computed only 
once. The cross sections need not to be approximated as in Fittino. I do not fit all MSSM 
parameters to the cross sections. I assume that all masses are known from LHC/ILC mea- 
surements. 
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Table A.I.: Vertex factors with parameters a, b, c, d, f, and g defined in Eqs. (|A.3|) , (|A.4|) , with 
e > 0. 





Vertex 


Factor 









Sl 


- <°' 




7~ 






7~ 




ie{pi + P2Y 


z- - 


< 


ieY {cPl + rfPs) 


z- - 







A. Radiative Neutralino Production 

A.l. Lagrangian and Couplings 

For radiative neutralino production 

e-(pi) +e+(p2) ^ fi{h)+x\{k2)+l{q), (A.l) 
the SUSY Lagrangian is given by [1] 

C = V2eaf,Pix°iSR + V2ebfePRf^ei + ^eZ^£^j>'[gPi+fPR]f^+K (A.2) 
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with the electron, selectron, neutralino and Z boson fields /g, Si^r, x^, and Z^,, respectively, and 
Pl,r = (l T 7^) /2. The couplings are 



cos 0„. 11' 

1 



2 sin 6; 



ig-(Ni2 + tan0,„Nii), 



(|Ni3|2-|Ni4p), / 



d 2 sin 6„, COS c/jy 

see the Feynman rules in Tab. lA.ll The Z-e-e couplings are 



1 (1 

sin 9u, cos 9j„ \ 2 



sin^0i„), d = —tanOzv 



(A.3) 



(A.4) 



A. 2. Amplitudes for Radiative Neutralino Production 

1 define the selectron and Z boson propagators as 

1 



m| - 2mlo - 2ki-k2 - ilzmz 

"■1 



(A.6) 



The tree-level amplitudes for right selectron exchange in the i-channel, see the diagrams 1-3 in 
Fig. 14. 1[ are 



Ml = 2ie^\af 

Ml = 2ie^\af 
M3 = 2ie^\af 



2pi-q 



u{ki)PRu{pi)] \v{p2)PLvik2)] {2pi - 2ki - q)-e* Ae^{pi,ki)Ae^{p2,k2) 



(A.7) 
(A.8) 

0- 
(A.9) 



The amplitudes for w-channel Sr exchange, see the diagrams 4-6 in Fig. 14.11 are 

Ae^(p2,fcl), 



M4 = -2ie'\a\^ 

Ms = 
Me = -2ie^\a\^ 



2ie^\a\^ 



u{ki)PR 2p^ q ^*u{pi)\ [v{p2)PLv{ki] 

u{k2)PRU{pi)\ [v{p2)r PLV{ki) Ae,{p,,k2), 



(A.10) 
(A.11) 



fi(fc2)PRw(pi)J \v{p2)PLv{ki)^ {2pi-2k2-q)-e* Ae^{pi,k2)Ae^{p2,ki). 

{A.12) 



In the photino limit, my amplitudes Mi-M(„ Eqs. (|A.7|| - (|A.12|) . agree with those given in Ref. [48]. 
The amplitudes for Z boson exchange, see the diagrams 7 and 8 in Fig. 14. 1[ are 



Mj = ie' 
Mb, = ie^ 



v{p2W{cPL + dPR) r»(pi)j \u{k,)j,,{gPL+fPR) v{k2)]Az{h,k2), 



2prq 



HP2)^* ^i (cPi + dPR) m(pi)1 \u{ki)jf, [gPi + /Pr) v{k2)]Az{ki,k2 

zp2-q -I L J 



(A.13) 

k2). 
{A.14) 
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Note that additional sign factors appear in the amplitudes M./i-M(, and Mj-Mg, compared to 
Aii-Ais- They stem from the reordering of fermionic operators in the Wick expansion of the S 
matrix. For radiative neutralino production e+e^ XiXil' such sign factors appear since the 
two external neutralinos are fermions0 For details see Refs. [38,58]. In addition an extra factor 
2 is obtained in the Wick expansion of the S matrix, since the Majorana field Xi contains both 
creation and armihilation operators. In my conventions I follow here closely Ref. [38]. Other 
authors take care of this factor by multiplying the Zx^Xi vertex by a factor 2 already [1]. For 
more details of this subtlety, see Ref. [1]. 

The amplitudes Aig — Aiu for left selectron exchange, see the diagrams 9-14 in Fig. 14.1 [ are 
obtained from the amplitudes by substituting 



a^b, Pl^ Pr, Pr ^ Pl, 
For Sl exchange in the t-channel they read 



= 2ie^\b\^ 
Mn = 2ie^\b\^ 



Mg 



u{h)PL 



2p2-q 



(A.15) 

(A.16) 
(A.17) 



m(/:i)Pi.m(pi)J [v{p2)PRv{k2)^ {2pi - Iki - q)-e* Ae^{pi,ki)AedViM ■ 

{A.18) 



The i/-channel exchange amplitudes are 

Mi2 = -2ie^\b\^ 



u{k2)PL ^l\ r<pi> 



2pi-q 



v{p2)PRv(ki) As^^{p2,ki), 



Mi3 



-2ie'\b\^ [ll{k2)PLU{pi)] \v{p2)^*^{^^PRV{k,) Ae,{pi,k2), 



2pi-q 



(A.19) 
(A.20) 



-2ie3|&|2 



u{k2)PLU{ 



v{p2)PRv{ki)^ {2pi-2k2-q)-e*Ae,ipi,ki)Ae,ip2,k2). 

(A.21) 



Our amplitudes Mi-Aiu agree with those given in Ref. [57,58], however there is an obvious 
misprint in amplitude T5 of Ref. [57]. In addition I have checked that the amplitudes Aii = 
e^Mf for z = 1, . . . , 14 fulfill the Ward identity q^ (E; M^) = 0, as done in Refs. [54, 57]. I find 
qjiiM'^ + M2 + A^s) = for f-charmel Cr exchange, q^iM^^ + M'^ + M'^) = for w-channel 
Sr exchange, qj, {My + -Mg) = for Z boson exchange, and analog relations for the ei exchange 
amplitudes. 



A. 3. Spin Formalism and Squared Matrix Elements 

I include the longitudinal beam polarisations of electron, P^-, and positron, Pg+, with — 1 < 
Pg± < +1 in their density matrices 



(A.22) 



^ Note that in Ref. [54] the relative sign in the amplitudes for Z boson ar\d t-channel Sr exchange is missing. 
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where cr^ is the third Pauli matrix and A_, A'_ and A+, are the helicity indices of electron and 
positron, respectively. The squared matrix elements are then obtained by 



'■<;■ 

where an internal sum over the helicities of the outgoing neutralinos, as well as a sum over the 
polarisations of the photon is included. Note that I suppress the electron and positron helicity 
indices of the amplitudes Aif^^* in the formulas (|A.7|) - l|A.2Hl . The product of the amplitudes 
in Eqs. (IA.24I) and (IA.25I) contains the projectors 

m(P,A_)m(p,A'_) = ^(^A_A'_ + 7'^_A'_)^ (^-27) 
vipA'+Mp.A^) = ^(VaV + 7'<aJ^- (A.28) 
The contraction with the density matrices of the electron and positron beams leads finally to 

pI v{p,A'+)v{p,K) = fi^PL + ^^PR^^ (A.30) 

A+,AV V ^ ^ / 

In the squared amplitudes, I include the electron and positron beam polarisations in the coeffi- 
cients 

Cr = (1 + P,-)(1-P,+ ), Cl = (1-P.-)(1 + P.+ ). (A.31) 

In the following I give the squared amplitudes Tij, as defined in Eqs. (|A.24|I and (|A.25|) . for cr 
and Z exchange. To obtain the corresponding squared amplitudes for ei exchange, one has to 
substitute 

b, d ^ c, f ^ g, Cr^ Cu A^^ Ag-j. . (A.32) 

There is no interference between diagrams with Sr and ei exchange, since I assume the high 
energy limit where ingoing particles are considered massless. 
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11 



T22 



T44 

T55 

T77 



(A.33) 
(A.34) 



T33 = 4e6c^|^|4^2^(p^^^^)^2^(p^^^^) p^.^^ p^.^^ _4m2^+8prfci +4(?-pi -4(?-/ci 



(A.35) 
(A.36) 



(A.37) 



4e^CR|fl|^ A?^(pi,/c2)A?^(p2,/ci)pi-fc2 P2-fci -4m2o +8pi-fc2 +4^j-pi -4(?-fc2 



{A.38) 



46'' 

'?-Pi 



4e6 
^■P2 



|Az(^i,^2)|' (Cr^V' + ClcV) n-h q-ki + {CR(f^ + Cic^f) Pi-ki q-h 



2^2 I r ^2r2\ 



Al 



(A.39) 



Az(^i,^2)r (CRd'f + ClcY) Pr^2 ^-^i + (CsdY + Cl^T) Pr^i '?-fc2 



-fg{Cic +CRd )m^o q-pi 



(A.40) 



T12 = -^e^CR\afAe^{pi,k^)Ae^{p2,k2)- 



13 



T 



14 



'■Pl i?-p2 

i?-fc2 Pr^i prP2 - pr^i '?-p2 Prfc2 + pr^i P2-fc2 'J-pi + PrP2 q-h P2-k2 
-q-pi p2-k2 pi-ki + pi-ki p2-k2 q-p2-2pi-p2 pvki p2-/c2 



(A.41) 



8e'CR\afAe,{Pi,h)Al{p2,k2)^^ 
m\q-pi +2(pi-fci f+ pi-ki q-pi -Ipi-ki q-ki 



m\q-p2 
-4e6CR|fl|4A,-,(p2,/ci)A,,(p2,/c2)- 



q-pi 



(A.42) 
(A.43) 
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^ Pi-Pi r 
4e^CR\a\^AeM'k2)AeM'k2)— 

q-pi q-p2 



q-pi - pi-pi + q-p2 



(A.44) 



4e^CR|fl|^Ae^(pi,fc2)Ae-^(p2,fcl)Ae-^(p2,fc2) 



q-pi 



-2pi-/c2 pi-pi - q-pi pi-pi + q-k2 pvpi - q-pi pi-k^ + q-p^ pi-k2 (A.45) 
4e^\a\^CRd-^A,^{p2,k2)mt{Az{ki,k2)}[2gp2-k2 q-ki + fm^^oq-p2] (A.46) 



-4e^CR\a\^d- 



-A,,(p2,/c2)9^e{Az(/ci,/:2)} 



■pi q-p2 

g(^-lpi-p2 p2-k2 pi-ki + p2-k2{q-ki pi-p2 - P2-ki q-pi + pi-ki q-p2) 
+ Pi-kiiq-pi p2-k2 + q-k2 pi-p2 - q-p2 Prfc2)) 



-fm^opi-p2{pvP2 - q-p2 - q-pi 



8e^CR\a\ 



4 Prh .2 



q-P2 



Al ipi,ki)^eAPl'^l) 



m\q-p2 + 2{p2-k2f+ ^2-^2 q-p2 -2p2-^2 q-k2 



m^^oPi-p2 

4e^CR|fl| Ae^{piM)Ag^{p2,ki)^^^-^[ q-pi - PvP2 + q-p2 



{A.47) 



(A.48) 
(A.49) 



4e^CR\a\''AsApiM)AeAPiM)- ' 



q-pi 



'^e^CR\a\^Ae^{pi,k2)Ae„{p2,ki)Ae„{pi,ki) 



2 



q-p2 



-lp2-ki pi-p2 - q-p2 PrP2 + q-ki prp2 - '?-p2 Pr^i + t?-?! Prki (A.51) 



Ae^CR\a\^d 



A,,(pi,fci)5He{Az(/ci,/c2)} 



(?-pi (?-p2 

g[lpi-p2 p2-k2 pi-ki + p2-k2{-q-ki prp2 + p2-fci (?-pi - prh (?-p2) 



+ pi-ki{-q-pi p2-k2-q-k2 prp2 + ^•p2 ^1-^2) 



{A.50) 



-fmloPvP2{pvP2 - q-p2 - q-pi) 



(A.52) 
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T: 



28 



T34 = 



— -^±-^Ag^{p^,ki)^i{Az{h,k2)}[2gpi-ki q-ki +/mJo'?-pi 
-4e^CR\a\^^^AeM'ki)AeM'ki)K~AP2,k2) 



lpi-p2 pi-ki + pi-p2 q-pi - pi-ki q-p2 + p2-ki q-pi - pi-p2 q-h 



(A.53) 



(A.54) 



T35 



2 

m 

-4e^CR\af^^Ae^{pi,ki)Ag^{pi,k2)Ae^{p2,k2) 
q-p2 

IprPi prk2 - pvP2 q-k2 + PrPi q-pi - P2-'<2 q-pi + Pvk2 q-pi 



(A.55) 



T36 = %e^CR\a\^Ae^{pl,kl)Ae^{pl,k2)^et.{pl,kl)^el.{P2,k2) 



T37 



T38 



T45 



m^o PrPi 



Ae^CR\a\^d 

q-pi 



-Ipi-ki -Iq-pi -2pi-k2 +2ki-k2 + q-k2 + q-ki 



A,-^ ( pi , fci ) ( P2, ^2 ) { Az (fci , /C2 ) } 



Igp2-k2[q-p\ pvki -2pi-kiq-ki +2{pi-kif + m^^oq-pi) 



(A.56) 



+/mJo(2prp2 pi-ki + prP2 q-pi - PvPi q-ki - pvki q-p2 + q-pi p2-ki){A.57) 



4e^CR\a\^d 
q-p2 



Ag^{pi,ki)As^{p2,k2)'nt{Az{ki,k2)} 



2gprki {2{p2-k2f+ prk2 q-pi -2p2-k2q-k2 +m^^oq-p2) 

h/m^o(2prp2 P2-^2 + Pi-pi q-pi - PvPi q-k2 + q-pi Pvk2 - q-pi p2-fc2)(A.58) 

{pi,k2)Ae„{p2,h) 

q-pi q-p2 

q-ki pi-k2 pi-p2 - Pvk2 q-pi pi-ki + prfc2 p2-ki q-pi + prp2 q-k2 Pi-ki 



q-pi p2-ki p2-k2 + prk2 p2-ki q-p2 -2prp2 Pr^2 P2-'^i 



(A.59) 



746 



Se'CR\at^^eAPl,k2)^Up2,k^) 
q-pi 



m\q-pi +2(prfc2) + Pi -^2 q-pi -2pvk2 q-k2 



(A.60) 
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4s ^ C I fl P c/ r 1 

^^A,-Jp2,fci)9^e{Az(fci,/c2)} 2gp2-h q-ki + fm\q-p2 (A.61) 

q-pi L ^1 J 

hpi-k2{-q-pi pi-ki + (?-p2 pi-fci - q-ki pi-pi)^ 



(A.62) 



Vi-ki 
q-p2 



A?Jpi,fc2)A,-^(p2,fci) 



p2-fci (?-p2 -2p2-ki q-ki +2(p2-fci) +m oq-pi 



(A.63) 



4e6CR|fl|2d 



A,,(pi,fc2)^e{Az(fci,/c2)} 



(J-P2 (?-pi 

g(2pi-p2 pi-k2 prki + pi-k2{-pvVi q-ki + pr^i (7-p2 - q-pi pi-ki) 
Vp2-ki (-prp2 '?-fc2 - Vi-ki q-p2 + (J-pi ^2-^2) 



-f^\oPrPi{prPi - q-p2 - q-pi 



4e^CR\a\^d 
q-p2 

Ae^CR\a\^d 

q-pi 



Aej,(pi,fc2)5^c{Az(fci,/c2)} [2gprfc2 q-ki +fm^^oq-pi 



A,-Jpi,fc2)AeJp2,fcl)^e{Az(/Ci,/C2)} 



{A.64) 
(A.65) 



2gp2-ki {pi-k2 q-pi -lq-k2 prk2 +2{pi-k2f + m^^^q-pi) 

h/w^o(2prfc2 prp2 + q-pi prp2 - q-k2 prpi - q-pi Pr^2 + q-pi p2-fc2)(A.66) 



q-p2 



- A,-Jpi, /:2) A,-Jp2, fci )^Ke{ Az (fci, /C2)} 



2gpi-k2 (2(p2-fci)^+ q-p2 pi-ki - 2p2-ki q-ki 



m 



q-Pi) 



-/mJo(2prp2 P2-fci + Pi-pi q-pi - Pi-h q-pi + pr/ci q-p2 - PrPi q-ki)(A.67) 



4e^ 



q-p2 q-pi 



— |Az(fci,/c2)P[(CR/rf' + CL/V)(2prp2 pr/ci p2-fc2 
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+ Pi-h {pi-k2 q-lpi - pi-p2 q-ki 
+ p2-k2 {p2-h q-pi - pvp2 q-h 
+ {CLgh^ + CRfd^)(2pvp2 Pvk2 pi-h 

+ P1-/C2 (pr^i q-p2 - Pi-P2 q-ki 
+ Prki {p2-k2 q-pi - prpi q-k2 

+2gf{CLC^ + CRd^)m^^opi-p2{pi-p2 - q-p2 - q-pi) 

I have calculated the squared amplitudes with FeynCalc [122]. When integrating the squared 
amplitude over the phase space, see Appendix the s-i-interference terms cancel the s-ii- 
interference terms due to a symmetry in these channels, caused by the Majorana properties 
of the neutralinos [55]. Note that in principle also terms proportional to e^xyiyk\p\p2'f3m{ls.z} 
would appear in the squared amplitudes T^, due to the inclusion of the Z width to regularise 
the pole of the propagator Az. However, since this is a higher order effect which is small far 
off the Z resonance, I neglect such terms. In addition they would vanish after performing a 
complete phase space integration. 



P2-k2 q-p2 
Pvki q-pi 



P2-ki q-p2 
Pvk2 q-pi 



(A.68) 
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B. Amplitudes for Radiative Neutrino 
Production 



For radiative neutrino production 

e-(pi) +e+(p2) ^ t/(fci) +17(/C2) +7('?). 
I define the W and Z boson propagators as 

1 

Aw{pi,kj) 



mlf + 2pi-kj 
1 



m\ — lki-k2 —iTzmz 

The tree-level amplitudes for W boson exchange, see the diagrams 1-3 in Fig. lB.lj are then 



Ml 



Ml = 



^wipiM) v{p2)YPLv{ki) u{ki)j^,Pi{^ - p^)^*u{pi) 



Aw{pi,ki) u{kl)J^'PLu{pl) v{p2)f{f^-i^)^^iPiv{k2] 



M^ = -ie^a^Awipi,ki)Aw{p2,k2) u{ki)j^PLu{pi) vip2)7^PLvik2) 

{{2ki - 2pi + q)^ig^ii + (pi -ki- 2q)figj,^ + (pi - k^ + q)c^g^^,) (e'') 
with the parameter 

1 

a = — . 

sm 

The amplitudes for Z boson exchange, see diagrams 4 and 5 in Fig. IB.ll are 

ieV 



M4 



with the parameters 
1 



Az(/Cl,fc2) \u{ki)YPLV{k2)\ [v{p2hv{cPL + dPR)i^ - f^)fu{pi) 

Az{k^M)\u{k{)7"PLo{k2^ \v{p2)r{f2-(^)lv{.cPi^^dPR)u{p{) 



sin Q-,0 cos 07„ V 2 



1 2 A 

- - sin 1 , d = - tan / 



sm v-,n cos 



(B.l) 

(B.2) 
(B.3) 

(B.4) 
(B.5) 

(B.6) 
(B.7) 

(B.8) 
(B.9) 



(B.IO) 



I have checked that the amplitudes Mi = ^itM^ for / = 1, ...,5 fulfill the Ward identity 
qpiiLiM'-) = 0. I find + Ml + Ml^) = for W exchange and q,r{Ml + X^') = 

for Z exchange. 
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e 

W+i 

^ L_ 



W+i 

e 



W+i 

e " ' — *—Ve 




< 



< 



diagr. 1 



diagr. 2 diagr. 3 diagr. 4 diagr. 5 

Figure B.I.: Contributing diagrams to e^e^ — > vv'y [77]. 



Table B.I.: Vertex factors with the parameters a, c, d, and / defined in Eq. (|B.7|I and (|B.10|I . 



Vertex 



Factor 











7 < 


-ie[(fci - + (/C2 - fc3)/J^,/a + (^3 - ^\)}igc,^\ 






w+-.-^ 


1 



I obtain the squared amplitudes T„ and T,, as defined in Eqs. (|A.24|) and (|A.25|) : 



Til = 



T22 — 



Aw(Pi'^i) Pr^2 q-ki 



q-p2 

T33 = e^CLfl''Afv(p2,^2)Afv(pi,/Ci) p2-fc2 Pl-fcl Pl-fcl + (p2-fcl (7pr/C2 -6(?-/C2) + 

prkiiq-h - q-pi)- q-ki ( prP2 + 2(?-p2 ) + ^1-^2 ( PrP2 + 6i?-p2 )) Pr^i + 
p2-fci q-ki prfc2 - 3q-ki prk2 prpi + q-ki q-k^ prp2 - p2-fci Vrki q-Tpi + 
(^■fci P2-/C2 q-fi +2p2-fci (?-fc2 '?-pi +2q-ki pi-/c2 (?-p2 + /ci-fc2 {-2q-ki pi-p2 - 

pi-fci (p2-fci - prP2 + q-p2)+ q-pi {5p2-ki +2pi-p2 + q-pi)) 



(B.ll) 
(B.12) 



'44 



'55 



'?-Pl 

q-p2 



\^z{ki,k2)\ {Cl^ p2-ki q-k2 + CrcI^ p2-k2 q-ki) 



\^z{ki,k2)\^{CiC^ pi-k2 q-ki +Qd^prfci q-k2) 



(B.13) 
(B.14) 

(B.15) 
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2p2-ki pi-k2 pi-p2 - q-h pi-fcs pi-pi - prh q-ki prpi - prh p^-ki q-pi + 



p2-ki p2-k2 q-pi + pi-ki pi-k2 q-p2 - Pi-ki pi-k2 q-p2 



(B.16) 



q-pi 



Aw{p2,k2)^w{pi,h) 



4pi-fci p2-ki pi-k2 - p2-ki q-ki pi-k2 - 3q-ki pi-p2 pi-k2 + 3prki q-p2 pi-k2 
3pi-ki p2-ki q-k2 + q-ki q-k2 pr/?2 + ^1-^2 Pi-ki q-pi - pi-ki p2-k2 q-pi + 
3p2-ki q-k2 q-pi + ki-k2 pi-p2 q-pi - Pvki q-k2 q-p2 

le^Cicfa^ 



q-pi 

e^Cicfa^ 



■Aw(p2,fc2)9^e{Az(fci,/c2)} Pi-ki q-k2 
■Aw{p2,k2)9\t{Az{ki,k2)} 



(B.17) 
(B.18) 



q-pi q-pi 

2p2-ki pi-k2 pi-p2 - q-ki pi-k2 pi-p2 - pi-ki q-k2 pi-p2 - 
Pi-h pvk2 q-pi + p2-h P2-k2 q-pi + Pr^i Pi-k2 q-pi - pi-ki pvk2 q-pi 



(B.19) 



q-p2 



Aw{Pi'h)Aw{p2,k2) 



-3pi-k2 p2-ki p2-ki +3q-ki pi-fcz p2-ki - pi-fci p2-ki p2-^2 + h-h PrPi Prh + 

2pvk2 prPi Pi-h - q-k2 prPi Prh -2pvk2 q-pi p2-ki + ^2-^2 q-pi Pi-ki - 
3pvk2 q-p2 pi-ki + pi-fci q-ki p2-k2 - ki-k2 q-ki prp2 - q-ki q-k2 prPi + 

q-ki p2-k2 q-pi +2pi-ki prk2 q-p2 +3q-ki prk2 q-p2 (B.20) 

q-pi q-pi 

2p2-ki pi-k2 pi-p2 - q-ki pi-k2 pi-p2 - Pi-ki q-k2 pi-p2 - Prki pi-k2 q-pi + 
p2-h Prk2 q-pi + Pi-ki pi-k2 q-p2 - Pi-ki pi-k2 q-p2 
2e^CLcfa^ 



q-pi 

e^Cicfa^ 

q-pi 



-Aw(pi,fci)9^e{Az(fci,fc2)} Pi-k2 q-ki 



Aw(pi,fci)Aw(p2,fc2)^e{Az(fci,^2)} 



(B.21) 
(B.22) 



4pi-/ci p2-ki pi-k2 - p2-ki q-ki pi-k2 -3q-ki prp2 pr^i +3pr/ci q-p2 pr^i 
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T35 



T45 



3pi-ki p2-h q-k2 + q-h q-k2 pi-pi + h-k2 p2-ki q-pi - pr/q p2-k2 q-pi 
3p2-ki q-k2 q-pi + krk2 prpi q-pi - Pr^i q-k2 q-pi 
e^Cicfa^ 



(B.23) 



?-P2 



-Aw{pi,ki)Aw{p2,k2)^e{Az{ki,k2)} 



-3pi-k2 p2-ki p2-ki +3q-ki pi-k2 prki - pi-h Prk2 pi-ki + fcr/c2 pi-p2 pi-h + 
2pi-fc2 prPi Pi-ki - q-k2 prPi Prki -2prk2 q-pi pi-ki + prk2 q-pi pi-ki - 
3pi-k2 q-p2 Pi-ki + pi-ki q-ki ^2-^2 - ^1-^2 PrPi q-ki - q-ki q-k2 pi-p2 + 
q-ki p2-k2 q-pi +2pi-ki pi-/c2 q-p2 +3q-ki pi-k2 q-p2 (B.24) 

3e^f 



q-pi q-pi 

.2 



\^z{kl,k2)Y 



Cic^{2p2-ki pi-k2 pi-p2 - q-ki pi-k2 prp2 - pi-ki q-k2 prp2 - 
Pi-ki prki q-pi + pi-h p2-^2 q-pi + Pvh pvk2 q-p2 - Pi-ki pvk2 q-p2) + 
CRif{2pi-ki prk2 prPi - q-h pi-k2 prPi - Pvh q-k2 prp2 - 
Pi-h Pi-ki q-pi + pi-ki p2-k2 q-pi + pi-h pi-k2 q-p2 - Pvki p2-k2 q-pi) 



(B.25) 



I have calculated the squared amplitudes with FeynCalc [122]. I neglect terms proportional to 
£K\}ivk\p\p2q"3m{Az} , see the discussion at the end of Appendix lAl 
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C. Amplitudes for Radiative Sneutrino 
Production 



For radiative sneutrino production 

I define the chargino and Z boson propagators as 

1 



ml + 2pi-kj 



m\ — 2m? — 2 fci ■ ^2 — iTzntz 



(C.l) 

(C.2) 
(C.3) 



The tree-level amplitudes for chargino x.^ exchange, see the contributing diagrams 1-3 in 
Fig. lC.ll are 



Ml 
Ml 

M3 



11 



}^0^A^,{p,,k,) [v{p2)nP2 - ^)Mh -fi- m^)Vi^u{n)\ (C.5) 



-ie^a^\Vn\^A^^{piM)\^{p2,k2) 



{C.6) 



with the parameter a defined in Eq. (|B.7|) . The 2x2 matrices U and V diagonalise the chargino 
mass matrix X [1] 



U*XV 1 = diag (ni^+, 



{C.7) 



The amplitudes for chargino x.2 exchange, see the contributing diagrams 4-6 in Fig. lC.ll are 

ie3fl2|y2i|2 



M^ 

Ms 
Me 



2q-pi ^2 

ie3fl2|y2i|2 
Iq-pi 



A^+(P2,/C2) [v{p2)PR{f2 - ^2 - m^^)Pi^{f^ - (jj^^Cpi)], (C.8) 
\i{Pi,h) [v{p2)r{f2 - 9)Mh - Fi - ^xi'^PMPi)]^ (C.9) 



-ie'a^\V2iW^ipi,ki)A^^{p2,k2) 

v{P2)Pr{P2 -^2- m^+)^*{Ki - - m^+)PMPi) 



(C.IO) 
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-f\/'\y\y\,'y 



diagr. 1 



xt 



diagr. 2 



xt 



xt 



xt' 



xt' 



-7 

-■V* 



diagr. 3 



xt' 



diagr. 4 



.■■■V 

■V* 




diagr. 5 



diagr. 6 



diagr. 7 



diagr 



Figure C.I.: Contributing diagrams to e^e — > vv*^ [77]. 



The amplitudes for Z boson exchange, see the diagrams 7 and 8 in Fig. IC.ll read 



Mi 



4q-p2 



Az(fcl,fc2) \v{p2)r{g - F2)(^l - ^2){cPl + dPR)u{p,] 



(C.ll) 



(C.12) 



with the parameters c, d, and / defined in Eq. dB.lOl ). I have checked that the ampHtudes Aii = 
e^M^l, z - 1, . . . , 8, fulfill the Ward identity q^iX^i -M f) = 0, as done in Ref. [81]. I find q^^ {M\ + 
Ml + Ml) = for exchange, q^t{Ml + M^, + M^) = for exchange, and q^^iM^ + 
M^) = for Z boson exchange. Our amplitudes for chargino and Z boson exchange agree 
with those given in Refs. [81, 82], and in the limit of vanishing chargino mixing with those of 
Ref. [52]. However, there are obvious misprints in the amplitudes M2 and M4 of Ref. [82], see 
their Eqs. (7) and (9), respectively, and in the amplitude T5 of Ref. [52], see their Eq. (F.3). 
I then obtain the squared amplitudes T„ and T,-, as defined in Eqs. (|A.24|) and (|A.25|I : 



Til = 



T22 



e^C^fl'^lViil^ 

— 2^.p/^ A^+(P2,/c2)(2p2-fe2 q-'ki -mlq-p2 



-A^+(pi,fci)(2pi-fci q-ki - mlq-pi 



2q-p2 

T33 = e^CLfl'^|yii|'^Aj+(pi,A:i)Aj+(p2,^2)[wJ+PrP2 +^m^^+prki p2-k2 - 
Imlpi-ki p2-ki +4ki-k2 pi-ki ^2-^2 -2m?pi-fc2 ^2-^2 +mj,pi-p2 



44 



— - — ' ^Ar,+ (p2,^2)(2p2-fc2 q-h - m^q-p2 

Zq-pi X2 



T55 



21 



lq-p2 



A^+(pi,/ci)(2pi-fci q-ki - mlq-pi 



Tee = e^CLfl'^|V2i|'^A^+(pi,/ci)A^+(p2,/c2) PrP2 +4ra^+prfci p2-fc2 - 
Imlpi-ki p2-ki +4/crfc2 pr^i P2-k2 -2m?pr/c2 p2-fc2 +mjpi-p2 



(C.13) 
(C.14) 

(C.15) 
(C.16) 

(C.17) 
(C.18) 
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Table C.I.: Vertex factors with parameters a, f defined in Eqs. (|B.7|) and (|B.10|) . and C the charge 
conjugation operator. 



Vertex 



Factor 



e,}i,T 



-ieaVjiPRC, transposed 




T77 



T12 



[p2-ki q-ki - p2-k2 q-h - p2-^i q-k2 + P2-^2 '?-fc2 -mlq-p2 + /cr^z i?-p2^ 



■|Az(fcl,/C2)|^ 



[pi-fci - pi-k2 q-ki - pi-ki q-k2 + pi-k2 q-k2 - m^q-pi + krk2 q-pi 



(C.19) 



(C.20) 



q-pi q-p2 



A^+(pi,/ci)A^+(p2,^2) 



'^1 



T- 



13 



[- ^1-^2 pr/?2 Pi-P2 + Prki pi-k2 pi-p2 + pi-ki ^2-^2 prP2 - 

q-ki p2-k2 prPi - Pvki q-k2 prPi + ^1-^2 q-pi PrPi + ^1-^2 q-p2 PrPi 

P2-ki pi-k2 q-pi + p2-ki p2-k2 q-pi + pr^i pr^2 q-p2 - Prki ^1-^2 '?-p2. 



(C.21) 



-A.+ (pi,fci)A^^+(p2,fc2) 



14 



[m^^+q-k2 pvp2 + m^^+q-pi p2-^2 -m^^+q-p2 prk2 -4pi-fci pi-fc2 P2-fc2 
4pi-fci P2-/C2 (?-fc2 +2mlpi-ki pi-p2 - Imlpi-ki q-p2] 



q-vi 



A^+(p2,fc2)A^+(p2,fc2)[2p2-fc2 q-k2 -mlq-p2 



X2 



(C.22) 
(C.23) 
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C. Amplitudes for Radiative Sneutrino Production 



= - \tipiM)\.{p.h) 



e'CL^\Vnl\V2i\^ 

q-vi q-vi 

[ - ki-k2 pi-p2 PrP2 + P2-ki pi-k2 pvP2 + Pi-ki p2-^2 Pi-P2 - q-ki prk2 pvP2 - 

Pi-ki q-k2 pi-p2 + ki-k2 q-pi pi-p2 + ki-k2 q-p2 Pi-p2 - P2-ki pi-k2 q-pi + 

p2-ki p2-k2 q-pi + prh pvk2 q-p2 - P2-ki pvk2 q-p2] (C.24) 

Tie = - ''^'"'l^^pf A,, (pi. fei ) A^. {p2, k2)\. {P2, k2) 

[m^^+q-k2 pvp2 +m^+(?-pi p2-^2 -m^^+q-p2 Pvk2 -4pi-fci pi-/c2 p2-/c2 + 

4pi-fci P2-/C2 'J-Zq +2mlpi-ki pi-p2 - 2m?pi-fci (^•p2] (C.25) 

= - "'^"2'^';p|^" V fe2)ine{ Az (fei, ^2) } 

p2-k2 - 2q-k2 p2-k2 + P2-fci q-k2 + mlq-p2 - fcr/c2 '?-p2] (C.26) 

e<^CL^cf\Vn\^ 
2q-pi q-p2 

[ - q-p2 Pi-k2 Prh + p2-fci prP2 Pi-k2 - Ip2-k2 pvP2 Pvk2 + q-k2 pvp2 prh - 
Pi-ki q-pi pvk2 + p2-k2 q-pi Pvk2 + prki q-p2 Prk2 - P2-ki q-p2 Prh + 
P2-h i?-p2 Prh + mlpi-p2 prP2 - h-h PrP2 PrP2 + Prh P2-h PrP2 - 
q-ki p2-h PrP2 - Prh q-k2 PrPi + P2-h i?-h PrPi - P2-h P2-h q-pi + 
p2-h P2-h q-pi -mlprp2 q-pi + h-h prP2 q-pi -mlpi-p2 q-p2 + 
ki-k2 pvP2 q-p2] (C.27) 

T23 = -^^^^:^^Aj4Pi.h)A,.(p2,h) 
q-p2 ^1 ' Ai ^' 

[np-^+q-ki prp2 -m^+p2-h i?-pi + "2^+ Prh (?-p2 -4pi-h P2-h P2-h + 

4pi-h (?-h P2-h +2m?p2-h PrP2 -2m?p2-h (?-pi] (C.28) 

T24 = ^A +(pi,h)A„+(p2,h 

q-pi q-pi ^ ' Hi 

[ - h-h PrP2 PrP2 + P2-h Prh PrP2 + Prh P2-h PrP2 - q-h p2-h PrP2 - 

Prh q-k2 prP2 + h-h q-pi PrP2 + h-h '?-p2 PrP2 - P2-h Prh q-pi + 

P2-h P2-h q-pi + Prh Prh '?-P2 - P2-h Prh '?-p2] (C.29) 

= ''^''''|^^^'j''^'''' A^+(pi.h)A^+(prh)[2prh q-ki -mlq-pi] (C.30) 

r26 = ^-^^ LA^+(pi,h)A^+(prh)A^+(p2,h) 

[wp-^+q-ki pi-p2 -m^^+q-pi p2-h + m^^+q-p2 Prh -4pi-h P2-h P2-h + 
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4pi-ki p2-k2 q-ki +2m?p2-^2 PrP2 -2m?p2-fc2 q-Tpi] (C-31) 



e''Ciahf\Vn 



|2 



[(?-p2 Pr^i Vvh +2p2-ki pi-p2 pvh - q-ki pi-p2 pvh - p2-^2 VvVi Pr^i + 
q-k2 pvp2 Pi-h - Pi-h q-pi prh - Pi-h q-pi prh - pr'ki q-pi prh - 
mlprp2 VrVi + h-h VvVi PrVi - Vi'h q-h PvPi - Pi'h Pi-h PrPi + 
q-ki p2-k2 pvp2 + Pi-ki p2-ki q-pi + p2-ki pi-k2 q-pi - prh pi-h q-pi + 
m?prp2 q-pi - h-h PrPi q-p\ + Pi-h q-pi Pi'h +miprP2 q-pi - 
h-h Pi-P2 q-pi] (C.32) 

T28 = '-^^^^0^^A^.{piM)^^{Az{h,k2)} 

[ - q-ki prk2 +2q-ki pyki - pi-ki q-k2 - mjq-pi + fci-/c2 q-pi] (C.33) 

[m\^-{q-k2 pi-p2 + P2-^2 q-pi - Pi-h q-pi) -4pvh Pi-hipi-h - q-k2) + 
lmlprki{pi-p2 - q-pi)] (C.34) 



e^Ci^a^\Vu?\V2i 
q-p2 



|2 



T35 = A^+{pi,ki)A^+{pi,'ki)A^+{pi,h] 



[m^^+{q-ki prp2 + prki q-p2 - prh q-pi ) -4pvh prhiprh - q-ki ) + 
2mlprk2{pi-p2 - q-pi)] (C.35) 
136 = 2e6Ci.«^|yn|'|V2i|X+(pi,/:i)A^+(p2,/c2)A^+(pi,;ci)A^+(p2,/:2) 
[2pi-ki P2-/C2 (?«^+ +2^1-^2 ) + Pi -^2 - 

2m?(pi-fci Prki + pi-k2 prh) +mjpi-p2] (C.36) 



[mK{q-ki prp2 - q-k2 prPi + P2-fci q-pi - Prh q-pi - Pvh q-pi + 
pi-k2 q-p2) -2pi-h Prh Prh -2pi-ki pi-ki prk2 +2pi-ki prk2 q-ki + 
4pi-ki prk2 prh +2pi-ki p2-ki q-k2 -4prfci p2-h q-k2 -2m?prp2 prh + 
2ki-k2 Prh PrP2 +2m?pi-h q-p2 -2h-h Prh '?-p2] (C.37) 

= _^!^^|!|pl£A^+(prh)A,+ (p2,h)^e{Az(h,h)} 

[m\+{q-ki prPi - (j-h PrP2 - P2-h q-pi + Prh q-pi + Prh q-pi - 
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C. Amplitudes for Radiative Sneutrino Production 

Pi-h q-p2)+2pi-ki p2-'<2 Pi-h - ^pvh prh pi'h +^Pr'ki pi-ki q-h + 
2p2-h Pi-h Pi-h -2pi-k2 pi-k2 q-ki -Ipi-ki p2-k2 q-k2 +2m?prp2 ^2-^2 
2ki-k2 p2-k2 prP2 - 2mlp2-k2 q-pi + 2krk2 prk2 q-pi] 

I 211 A + A + p2,?C2 

q-pi q-p2 ^ ^ Az 



T45 



746 



T57 



(C.38) 



T47 = 



q-pi + 



[- fci-/c2 pvp2 PvP2 + Pvk2 P2-ki prPi + Pvh P2-k2 PvP2 - q-ki prk2 PvP2 
Prh q-h PrPi + h-h q-pi PrPi + ^1-^2 q-p2 PvP2 - P2-ki prh q- 
p2-ki p2-k2 q-pi + pi-ki pi-k2 q-p2 - P2-ki pi-k2 q-p2] 

2^ , 

q-pi ^2 ' X2 ^' 

[m^^+q-k2 pi-p2 +m^^+q-pi ^2-^2 -m^^+q-p2 ^1-^2 -^pi-h pi-k2 p2-k2 + 
4pi-fci p2-k2 q-k2 +2mlpi-ki prp2 -2m?pi-/ci q-p2] 

_ £!^^|M^ A^. {p2, k2)^^{^Z {h, k2)} 
[q-ki p2-k2 - 2q-k2 ^2-^2 + P2-ki q-k2 +mlq-p2 - ki-k2 q-p2] 

2q-pi q-p2 

[ - q-p2 Pi-k2 pi-k2 + p2-ki pi-p2 pi-k2 -2p2-k2 prp2 Pi-k2 + q-k2 prp2 Pi-k2 
Pi-h q-pi Pi-h + P2-k2 q-pi Pvk2 + prh q-p2 Prk2 - P2-ki q-p2 Pr^2 + 
k2 +mlpi-p2 pi-p2 - krk2 pvP2 PvP2 + Prh prk2 prP2 - 
Pi - Prh q-h PrPi + P2-^2 q-k2 PvP2 - P2-^2 ^2-^2 q-] 



(C.39) 



(C.40) 



(C.41) 



P2-k2 q-p2 Pi 
q-ki p2-k2 pi 
p2-ki p2-k2 q 
ki-k2 pvP2 q 



■pi + 



Pi 
Pi 



\prPi q-p\ + h-h PrP2 q-p\ -m\pYP2 q-p2 



(C.42) 



^ g'CLfl'|^2l|\ 2 , I, w , 

= \^{pi,h)A^^{p2,k2) 



q-p2 

[m^^+q-ki pi-p2 -m^^+p2-h q-pi Pi'h q-p2 -4pi-h P2-ki p2-^2 + 

4pi-^i q-ki p2-k2 +2mlp2-k2 pi-p2 -2m?p2-/c2 q-pi] (C.43) 

2q-pi q-p2 

[q-p2 Pi-ki pi-h +2p2-h Pi-P2 Pi-ki - q-h PrP2 Prh - ^2-^2 PrP2 Prki + 
q-k2 pi-p2 pi-h - Pi-k\ q-pi pvh - Pi-ki q-pi Pvh - Prh q-pi Pvki - 
mlprPi PrP2 + h-h Prp2 Pvp2 - p2-h q-h PrP2 - p2-h prh PrP2 + 
q-ki p2-k2 pi-p2 + P2-fci P2-h q-pi + P2-ki pi-k2 q-pi - p2-ki p2-k2 q-pi + 
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tnlpi-p2 q-pi - h-k2 pi-p2 q-pi + pi-h q-p2 pr^i +?«^PrP2 q-pi - 

h-'ki pvpi q-pi] {CM) 

[2q-ki pi-ki - q-ki pi-k2 - pvh q-k2 -rn^q-pi + ki-k2 q-pi] (C.45) 



[np-^+{q-ki pi-p2 - q-k2 prp2 + prki q-pi - p2-k2 q-pi - prki q-p2 + 
pi-k2 q-p2) - 2pi-ki Prki pi-k2 -2pi-ki pi-ki pi-^i +2prki ^2-^2 q-h + 
4prfci p2-'<2 Pr'<2 + 2pi-ki prki q-k2 -4pi-/:i p2-^2 q-^i -2m?prp2 Vvh + 
2^1-^2 prh VrPi +2m?pi-fci ^j-p2 - 2krk2 prh q-pi] (C.46) 

prp2 - '?-fc2 PrP2 - P2-fci '?-pi + ^2-^2 q-pi + Pr^i ^-^2 - 
prk2 q-p2) + 2pi-ki prk2 p2-k2 -4pr/ci p2-fci p2-fc2 +^prh pi-h q-ki + 
2p2-h Pi-h Pi-h -2pi-k2 p2-k2 q-ki -2pi-ki p2-/c2 q-'k2 +2m?pi-p2 ^2-^2 - 
2ki-k2 p2-k2 pi-pi - 2mlp2-k2 q-pi + 2ki-k2 p2-k2 q-pi] (C.47) 

Aq-piq-p2 

[pi-ki (pi-fci (?-p2 +2p2-/:i prp2 - q-h pvp2 -2p2-k2 prPi + q-'k2 PrPi - 
Pi-h q-pi + Pi-h q-pi - Pi-ki q-p2 - 2prk2 q-p2 + Prh q-pi) + 

pi-p2 ( -2m?prp2 + 2ki-k2 pi-p2 - pi-h q-h - 2p2-ki pi-k2 + q-ki pi-k2 + 
q-ki p2-k2 +2pi-/c2 p2-fc2 + P2-fci q-'k2 - Pi-k2 q-k2 - pi-h q-k2) + 

q-pi [pi-h Pi-h + P2-fc2 Pi-h + Pi-h Pi-h -2p2-ki ^2-^2 - 
pi-k2 p2-k2 +2m?pi-p2 -2ki-k2 prPi) + 

q-piipi-h Pr'ki + P2-fci Prh - Pvh Pi-h +2m?pi-p2 - 2krk2 PrPi)] (C.48) 

Formulas for the squared amplitudes for radiative sneutrino production can also be found 
in Refs. [81, 82] for longitudinal and transverse beam polarisations. Here, I give however 
my calculated amplitudes for completeness. I have calculated the squared amplitudes with 
FeynCalc [122]. I neglect terms proportional to eKAjii/^iPi P2^"^Tri{Az}/ see the discussion at the 
end of Appendix[A| 
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D. Definition of the Differential Cross Section 
and Phase Space 



I present some details of the phase space calculation for radiative neutralino production 

e-(pi) +e+(p2) - xl{ki)+xUk2) + 7{q). (D.l) 
The differential cross section for (|D.1|) is given by [123] 



d(7 



2 2s 



(D.2) 



where p f and E f denote the final three-momenta and the final energies of the neutralinos and 
the photon. The squared matrix element | ^ is given in Appendix [Al 

I parametrise the four-momenta in the center-of-mass (cms) system of the incoming particles, 
which 1 call the laboratory (lab) system. The beam momenta are then parametrised as 



-(Vi, 0, 0, -^s) 



(D.3) 



For the outgoing neutralinos and the photon I consider in a first step the local center-of-mass 
system of the two neutralinos. The photon shall escape along this Xj,-ax\s. I start with gen- 
eral momentum-vectors for the two neutralinos, boost them along their Xs-axis and rotate them 
around the Xi-axis to reach the lab system. Note that the three-momenta of the outgoing par- 
ticles lie in a plane whose normal vector is inclined by an angle 6 towards the beam axis. I 
parametrise the neutralino momenta in their cms frame [48] 



k* sin 6* cos cp' 
k* sin0* sincp'* 



rcosr / 

with the local cms energy s* of the two neutralinos 



( r.- 

r cos./)* 



V 



—k* siny^ cos(^' 
—k* sin0* sin^^ 
-k* cose* 



(D.4) 



J 



s* = {ki + k2f 



2m\ 



■2fci-fc2, 



(D.5) 



the polar angle 0* , the azimuthal angle (p* and the absolute value of the neutralino three- 
momenta k* in their cms frame. These momenta are boosted to the lab system with the Lorentz 
transformation 



(l 

10 

10 

V7iS 7/ 



(D.6) 
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with 7 = 



and ^ 



|ki+k2| 

I cms beam 



the boost velocity from the cms to the lab system 



(D.7) 



Boosting the momenta fc^ and k\, see Eq. (ID.4I) , at first with the Lorentz transformation Eq. (|D.6I > 
and then rotating with Q yields the neutralino and photon momenta in the lab system [48] 



7E* + fy^k* COS 0* 

k* sin 0* cos ^* 
/c* sin 0* sin ^* cos + ( j67£* + 7/c* cos 0* ) sin 
V -A;* sin 0* sin ^* sin + (,67E* + 7^* cos Q* ) cos 0/ 



^2 



with 



E* 
/37 

The differential cross section for e+e^ - 

1 



7E* - ^7^* cos Q* 
—k* sin0* coscp* 
-k* sin 0* sin ^* cos 9 + {^jE* - jk* cos 0* ) sin 6 
V k* sin 0* sin ^* sin 6 + (,67E* - 7A;* cos 0* ) cos / 







Ai 



2 ' 

s - s* 



2Vss* 
A'l/tiT riow reads [48] 



d(7 



40967r4s 

where the integration variables run over 

< (p* < In, 

-1 < COS0* < 1, 

4ra2o < s* < (l-x)s, 

-0.99 < COS0 < 0.99. 



1 



4m2n 



|7W|2dcos0dcosr d(^* ds* 



(D.8) 



(D.9) 



(D.IO) 



(D.ll) 
(D.12) 
(D.13) 



(D.14) 



(D.15) 
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E. How to calculate helicity amplitudes for 
longitudinal polarisation states 



E.l. Introduction 

Bouchiat and Michel presented [124] formulae to perform helicity spin sums for Dirac fermions 
and antifermions. Haber collected in [125] mathematical tools to deal with such sums and 
presented example calculations. Choi et al. [126] extended these formulae to spin 1 and spin 
I fields. In this paper I present a proof of the Bouchiat-Michel-formulae and I extend these 
formulae to Majorana particles. All formulae are written in a covariant manner. 

E.2. Spinor calculus 

The Dirac spinor s u{p,A) and v{p,A) obey the Dirac equation: 

{f — m)u = 0, {f + m)v = 0. (E.l) 

The charge-conjugation-operator C converts the spinor ti with positive energy into the spinor v 
with negative energy and vice versa: 

u = Cv^, V = CU^. (E.2) 

There are two solutions of the Dirac equations for a given 4-momentum ip, so there exists 
another good quantum number to label these states. This is the helicity A. The helicity operator 
A = Ep commutes with the Dirac operator ± m), so the eigenvalues A = it^ of A are good 
quantum numbers. 

The spinors u and v are normalized to ( [127]) 

u{p,A)u{p,A') = 2m5xA', (E.3) 

v{p,A)v{p,A') = -2m^AA'. (E.4) 

When calculating cross sections or decay widths of fermions and antifermions, I sum over all 
spins states and average over initial spins, using the completeness relation [127]: 

Y^u{p,A)u{p,A) = f + m, (E.5) 

A 

Y^v{p,A)v{p,A) = f-m. (E.6) 

A 

When describing spin-polarized fermion ensemble one introduces spin vectors. The longitu- 
dinal spin vector for a particle with mass m is defined by 

s is normalized to s ■ s = — 1, and is orthogonal to the momentum vector p: s ■ p = 0. 

I have to do distinguish between massive and massless particles. The spinvector for massless 
particles is obtained in the limit m 0. 
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E.3. The Bouchiat-Michel-Formula 



E.2.1. The massive case 

The operator 7^^ commutes with -p, so both operators can be diagonalized simultaneously. 
Their eigenvectors are known, and the eigenvalues are obtained by [128] : 

= 7^(^s- p-is^pi,cr'"'^i/(p,A) 
= -i7^ (cr^'sopj + a'°SjPo)u{p, A) 



= 2mEy-pj-M(p, A) 
= lXmu{p,\) 

= l\pu{p,\). (E.8) 

I have used the relation 7^7*^7^ = i'y0^i^2^3^oy _ _j^i^2^3y _ {^.^^ryk^i realize that 7^^ 
is the helicity operator. 

From the above calculation it follows immediately: 

7^" = 2Am. (E.9) 

E.2.2. The massless case 

The Dirac equation for a massless spin ^-fermion is: 

fu{p,k)=Q (E.IO) 

and multiply eq. dElOl) with 7^7° [128]: 

= ^S''Mp,k) = h^PQ-tp)u{p,X); (E.ll) 
tpu{p,K)=^hi{p,\) (E.12) 

with p° = \p\. The chirality-operator 7^ commutes with the helicity-operator, hence they have 
common eigenvectors, a similar equation holds for v spinors: 

-f^u{p,A) = ±m(p,A) = 2Am(p,A), (E.13) 
7^t;(p,A) = T^{p,\) = -2Av{p,\). (E.14) 

Since (7^)^ = 1, Trace(7^) = 0, the eigenvalues of the chirality-operator are ±1. 

E.3. The Bouchiat-Michel-Formula 

The Bouchiat-Michel-formulae (BMF) tell us how to contract spinors with different polarisa- 
tions. The BMF is interesting if the initial state fermions have the same mass, but also in using 
density matrix techniques [125]. 
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E. Helicity amplitudes 



E.3.1. Spin vectors 

I enlarge the set s, p defined in sec. IE .21 with two other four- vectors and to a orthonormal 
basis in space-time: 

p-s" = (E.15) 

s"-s^ = -S"^ (E.16) 

m 

with fl = 1 . . . 3, = s. 

The spinors u{p, A), v{p, A) satisfy: 

^'ru{p,A')=<4MP'^) (E.18) 

j'fv{p,A')=al,^v{pA). (E.19) 

The proof of this equation can be found in/follows [129]: Define the Pauli-Lubanski-vector as 

Wj, = ^e^^afiyM'^P'^ (E.20) 

= (fp, Et + Kxp) , (E.21) 



where M'*^, X' = M^', P'^ are the generators of the Poincare-group. In the rest frame the spin 
vectors take a simple form: s' = (O, e,) where e, is the ith unit vector. I build the scalar operator 
W • s\ which has the eigenvalues Am and evaluate it in the rest frame: 

W ■ s' = --mtci = --ma' (E.22) 

2 2 ^ ' 

Then, I apply W • s' on a spinor n{p,A): 

W-s'u{pA) = -^mcr{,^u{p,\'). (E.23) 

On the other hand W ■ s' = \ \f ,f\ 7^, so 

W-s'm(p,A) = -'^-{f,f\^\i{p,X) 
1 

= --^l^{ff-ff)u{p,A) 
= -^mj^^'u{p,A) 



2 

Now eq. (|E.18|) follows. The proof for the spinor v is similar. 

E.3.2. BMP for massive Dirac fermions 

Now I have all ingredients to formulate and prove the BMP: 



lm4,^u(p,A'). (E.24) 



u{p,A')u{p,A) = l[<5^,;, + 75^Vl,J(p + m) (E.25) 
v{p,A')v{p,A) = l\8^^,^^^fa\^\{f-m). (E.26) 
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E.3. The Bouchiat-Michel-Formula 



The sum must have the form 

u{p,\')u{p,X) = A5^,j, + B'al,^. (E.27) 
To determine the unknown coefficients A, B" we multiply both sides with S'^ and with o'^''^: 

^u{p,\')u{p,A)S^'^ = \[f^rn)=A, (E.28) 
1m(P,A')S(p,A)K)^'^ = U^f{^^m)=B\ 

(E.29) 

A similar proof holds for eq. ( IE.26I ) 

E.3.3. BMF for massless Dirac fermions 

To perform the limit m ^ I use (|E.12| |, and get 

u{p,\')u{p,\) = 

\ (^A'A + tM'a + + 7'^'4a) f' (E-30) 
v{p,k')v{p,k) = 

\ ('Jaa' + tMa' + tW],, + 7'^V2^,) f. (E.31) 

I shall make some remarks to the spin vectors for particles with E ^ ra. I start with eq. (|E.7|) 
and expand E = |p| y^l + m^/lpp ?a \p\ ^1 + + |^))- ^^^^ expansion preserves the 
normalization ■ = —1. 

E.3.4. Majorana-Fermions 

In supersymmetric (SUSY) field theories there appear Majorana fermions, for example the neu- 
tralinos, which are the SUSY partners of the neutral weak gauge and Higgs bosons. In Feynman 
diagrams with Majorana fermions I find often clashing arrows. So it is useful to have formulae 
to handle this case. I use the relations (|E.2|) , and I get 



uipA'Vip,^) 


= ^[^a'a + 7'«a](|^ + 








(E.32) 


v^ipA'MpA) 


= ^C-1[^a'a + 7'«'a]i 


[p + m], 


(E.33) 


v{p,A')u'{p,\) 


= ^[^a'a+7'«'a](^- 




(E.34) 


u^{p,A')v{pA) 


= ^C-1[Ja'a + 7'«'a]i 


[f - m). 


(E.35) 
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E. Helicity amplitudes 



E.4. Calculation of the density matrix 

The methods described above can be used to compute squared matrix elements in the helicity 
mechanism described in [125]. The only change in the usual mechanism is the completeness 
relation. I consider longitudinal polarized electrons. I treat the electrons as massless. I call this 
matrix the reaction matrix Raa'- 

Raa' = A)m(p, A') = ^(^^AA' + 7^(dA')r (E-36) 
The helicity indices are contracted with the beam matrix Baa'- 



Now I calculate: 



Baa'R'^' = ^(^'"' + PV3^^')1(^aa' + 7Ma> 



^^Pr + ^^Pl ) p (E.38) 



I can use eq (|E.38|) instead of the usual completeness relations Eqs (|E.5|l , (|E.6|) to do the spin 
sums. 
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